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I  INTRODUCTION 

This  is  the  final  report  of  a  project  devoted  to  the  theoretical 
study  of  molecular  state-to-state  collision  phenomena.  The  present 
document  describes  the  work  accomplished  during  the  fifth  year  of  the 
program  (1  May  1979  to  30  April  1980).  Interim  scientific  reports 
dated  June  22,  1979,  July  1,  1978,  June  24,  1977,  and  June  28,  1976 
have  covered  the  work  of  preceding  years . 

The  purpose  of  the  project  has  been  to  develop  the  theoretical 
tools  necessary  to  treat  a  variety  of  collision  processes,  and  to  per¬ 
form  calculations  of  transition  probabilities  for  selected  systems. 

The  problems  considered  have  involved  the  interactions  of  three  particles 
under  the  influence  of  potentials  for  two  or  more  electronic  states. 

For  this  type  of  process,  currently  available  experimental  techniques 
allow  the  determination  of  specific  states  of  species  before  and  after 
interaction.  Our  general  method  of  approach  has  been  to  pursue  two 
types  of  calculations:  model  calculations  designed  to  isolate  and 
probe  particular  reaction  mechanisms,  and  detailed  calculations  on 
specific  systems  that  have  been  studied  experimentally  in  our  labora¬ 
tory  or  elsewhere.  This  dual  approach  has  enabled  us  better  to  under¬ 
stand  a  wide  range  of  scattering  phenomena,  and  has  provided  a 
quantitative  test  of  the  accuracy  of  the  methods  we  have  developed. 

During  the  past  year  we  have  worked  on  the  following  types  of 
problems:  theoretical  studies  of  collisional  angular  momentum  mixing 
of  highly  excited,  Rydberg  atoms;  calculations  of  potential  surfaces, 
matrix  elements,  and  photodissociation  cross  sections  for  small 
molecules;  and  the  development  of  new  methods  to  treat  the  nonadiabatic 
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II  RESEARCH  ACCOMPLISHMENTS 


Theoretical  Studies  of  Rydberg  Atoms 

Studies  investigating  the  properties  of  highly  excited  Rydberg 
atoms  have  been  under  way  at  SRI  for  several  years.  The  experiments 
of  Gallagher,  Edelstein,  and  Hill^  first  identified  the  "angular 
momentum  mixing"  process: 

Na*(nd)  +  X  -  Na*(ni' )  +  X  (l'  St  3)  (1) 

in  1975.  The  cross  sections  for  this  process  are  very  large  and  show 

a  maximum  when  plotted  as  a  function  of  principal  quantum  number  n. 

2-5 

Over  a  period  of  years,  studies  supported  by  this  contract  have  con¬ 
tributed  to  increased  understanding  of  Rydberg  collision  processes. 

The  work  completed  in  the  past  year  has  furthered  this  understanding 
by  consolidating  the  results  of  several  previous  calculations  in  terms 
of  general  scaling  formulas.  These  formulas  lead  to  an  appealing 
physical  interpretation  of  the  process  of  angular  momentum  mixing,  and 
can  be  quickly  and  easily  applied  to  a  variety  of  systems. 

During  the  past  year,  in  work  supported  by  another  Air  Force 
contract,  the  research  group  of  T.  F.  Gallagher  has  studied  collisions 
of  Na  Rydberg  atoms  with  molecules  such  as  CH^  and  C^Hg  (methane  and 
propane).6  The  research  sought  to  determine  the  role  of  the  internal 
structure  of  the  collision  partner.  The  experimental  data  obtained 
was  qualitatively  very  similar  to  data  obtained  for  much  simpler 
collision  partners,  such  as  the  rare  gases  He,  Ne,  and  Ar,*’  and  also 
for  These  results  suggested  that  the  larger  and  more  complicated 

collision  partners  could  be  treated  by  the  same  theoretical  model^’^ 
that  was  successfully  used  to  analyze  collisions  with  structureless 
particles.  The  only  Information  required  for  this  model  is  data  for  low 
energy  electron  scattering  from  CH^  and  C3Hg,  which  is  available. 


Because  of  the  apparent  wide  range  of  validity  of  the  theoretical 
model  we  have  used,  and  because  of  the  small  number  of  parameters  on 
which  the  results  depend,  we  sought  to  express  the  results  in  terms  of 
physical  scaling  laws.  We  have  found  that  the  results  of  many  coupled- 
channel  and  Born  calculations  can  be  expressed  in  terms  of  numerically 
determined  functions  of  reduced  parameters.  The  general  result  is 

a£-mix(n)  "  ™  \  8(0)  f(v)  (2) 

where  a.  (n)  is  the  experimentally  measured  cross  section  for  process 
j2-mix 

(1),  summed  over  all  final  states  t!  ;  3  is  a  quantity  that  measures  the 
strength  of  the  interaction  of  the  collision  partner  with  the  diffuse 
electron  cloud;  and  y  is  a  quantity  related  to  the  energy  gap  between 
the  nd  and  higher  nl  states  of  Na.  Equation  (2)  thus  leads  to  an 
interpretation  of  the  scattering  process  in  terms  of  reduced  parameters 
that  have  a  clear  physical  meaning.  The  formula  has  been  applied  to 
collisions  of  Na(nd)  with  He,  Ne,  Ar,  N^,  CH^,  to  Rb(nf)  +  He,  Ar, 

Xe;  and  to  Xe(nf)  +  CC^.  The  results  are  generally  accurate  to  about 
a  factor  of  two.  This  degree  of  accuracy  is  very  useful  because  the 
cross  sections  vary  by  an  order  of  magnitude  or  more  in  the  range  of  n 
considered.  This  work  has  been  submitted  to  Physical  Review  A.  and  a 
preprint  is  included  as  Appendix  A.  This  preprint  contains  a  more 
complete  discussion  of  the  scaling  formula,  as  well  as  graphs  of  the 
functions  f(y)  and  g(3)  of  Eq.  (2). 

The  most  important  theoretical  result  of  our  work  on  Rydberg  atoms 
has  been  to  achieve  a  better  understanding  of  the  role  of  the  inter¬ 
action  of  the  excited  electron  with  the  collision  partner.  It  had 
generally  been  argued  that  the  interaction  of  a  Rydberg  atoms  with  a 
collision  partner  should  be  determined  by  the  way  the  excited  electron, 
in  its  large  orbit,  scattered  from  the  collision  partner.  The  quantita¬ 
tive  calculations  we  have  done  lend  support  to  this  physical  picture, 
and  explain  the  mechanism  in  much  greater  detail. 
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Theoretical  Calculations  of  Small  Molecules;  CH 


The  formation  and  destruction  of  CH+  in  diffuse  interstellar  clouds 
is  an  interesting  problem  in  interstellar  chemistry.  Better  knowledge 
of  the  excited  state  potential  curves  as  well  as  of  photodissociation 
cross  sections  is  necessary  to  understand  this  astrophysical  problem. 

The  CH+  system  is  small  enough  so  that  calculations  that  include  Rydberg 
states  on  the  same  footing  as  valence  states  could  be  undertaken.  In 
collaboration  with  Dr.  Bowen  Liu,  IBM  Research  Laboratory,  San  Jose,  CA, 
and  Dr.  Kate  Kirby,  Harvard -Smithsonian  Center  for  Astrophysics,  Cambridge, 
MA,  we  performed  ab  initio  calculations  of  the  valence  and  low-lying 
Rydberg  states  of  and  symmetry  of  CH+.  The  specific  objective  of 
this  work  was  to  explore  the  repulsive  excited  states  of  CH+  as  well  as 
to  look  for  bound  excited  states  that  are  dipole-connected  to  the  ground 
X^T+  state  to  investigate  whether  CH+  is  a  source  for  unidentified 
interstellar  molecular  absorption  lines.  Only  the  first  root  of  each 
symmetry  is  bound,  although  the  21£f  state  is  found  to  be  quasi-bound 
in  this  calculation.  Transition  moments  and  photodissociation  cross 
sections  from  the  X^2+  state  to  a  number  of  excited  and  states 
were  calculated;  it  was  found  that  the  3*Z+  and  21]!  states  have  signifi¬ 
cant  cross  sections  at  wavelengths  longer  than  Lyman-cr.  From  these 
cross  sections  and  astrophysical  information,  photodissociation  rates 
in  interstellar  clouds  and  comets  can  be  obtained. 

Fhotodissociation  cross  sections  and  their  astrophysical  implica¬ 
tions  are  discussed  in  a  paper  accepted  by  the  Astrophysical  Journal.^ 

The  Ah  initio  calculations,  in  which  care  was  taken  to  test  the  validity 
of  the  calculational  procedure,  are  described  in  a  paper  that  has  been 
accepted  by  the  Journal  of  Chemical  Physics.  A  preprint  of  the  latter 
paper  is  attached  as  Appendix  B. 
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We  have  also  recently  completed  a  calculation  of  the  £  state  of 

3  ^ 

CH  which  correlates  asymptotically  to  C(  P)  +  H  .  This  system  is  being 
studied  experimentally  in  the  Molecular  Physics  Laboratory  by  the 

technique  of  fast  beam  photofragment  spectroscopy.  Knowledge  of  the 

3  . 

£  potential  curve  may  aid  in  interpreting  the  experimental  data.  This 
particular  state  of  CH+  poses  an  interesting  theoretical  challenge  as 
well.  While  there  is  no  electron  on  the  H+,  asymptotically  there  is  H 
character  in  the  molecular  orbitals  in  the  bonding  region.  Thus,  de¬ 
signing  a  calculation  to  determine  properly  the  molecular  orbitals 
requires  some  care.  We  have  carried  out  a  series  of  calculations 
including  different  amounts  of  electronic  correlation  and  have  compared 
the  results  in  terms  of  the  spectroscopic  parameters  describing  the 
minimum  and  maximum  of  the  potential  curve.  The  excellent  agreement 

between  the  different  calculations  confirms  that  we  have  realistically 
3  * 

described  the  CH  £  potential  curve. 

Theory  of  Nonadiabatic  Collisions 

During  the  past  year  we  made  substantial  progress  in  the  study  of 
atom-molecule  collisions  that  involve  two  or  more  electronic  potential 
energy  surfaces.  Although  the  approaches  we  developed  are  applicable 
to  many  different  systems,  we  generally  focused  on  systems  that  were 
actively  being  investigated  in  the  experimental  program  at  SRI. 

We  have  used  three  types  of  theoretical  approaches.  The  first  is 

a  quasi-classical  trajectory  model  that  invokes  the  Landau-Zener 

formula  to  calculate  a  "surface-hopping"  probability  at  each  surface 

Intersection.  Such  a  model  has  yielded  reasonable  results  for  many 
8-11 

systems,  but  its  quantitative  accuracy  may  be  limited  in  certain 

1  3 

situations,  such  as  the  reaction  of  He  (2  *  S)  with  0^  now  being  studied 
at  SRI  —  a  reaction  that  involves  three  surfaces.  Therefore,  we  have 


7 


also  pursued  two  different  quantum  mechanical  approaches.  One  of  these 

is  based  on  applying  time-dependent  perturbation  theory  in  the  impact 

parameter  formulation.  We  assume  the  projectile  follows  a  straight 

line,  constant  velocity  classical  trajectory  of  a  given  Impact  parameter. 

The  potential  seen  by  the  target  diatomic  molecule  is  then  explicitly 

time-dependent.  The  other  quantum  mechanical  approach,  which  is  time- 

independent,  is  a  generalization  of  the  quantum  mechanical  Infinite- 

12  13 

Order-Sudden  (IOS)  approximation.  ’  This  method  has  been  extensively 

12  14 

applied  to  rotational  excitation.  *  In  this  approximation,  the 
solution  of  the  coupled  equations  is  greatly  simplified,  when  one 
assumes  the  rotational  levels  for  a  given  vibrational  state  are 
degenerate,  and  solves  coupled  channel  equations  for  only  the  number  of 
vibrational  states  of  interest.  We  have  found  that  similar  expressions 
can  be  derived  when  several  electronic  states  are  coupled.  They  can  be 
made  extremely  simple  if  one  assumes  that  both  the  rotational  and 
vibrational  levels  are  degenerate.  Although  this  approximation  must  be 
tested  further,  it  looks  promising  at  the  experimental  energies  of 
interest.  The  following  subsections  discuss  our  results  with  each  of 
these  three  approaches . 

(a)  Quasi-Classical  Surface  Hopping  Model 

The  work  done  in  this  area  was  an  outgrowth  of  a  collaboration 
with  Dr.  Aart  Kleyn  of  the  FOM  institute  in  Amsterdam.  Kleyn  visited 
SRI  for  two  months  in  1979  and  implemented  seme  of  the  computer  codes 
used  in  Amsterdam.  We  have  done  additional  work  to  modify  and  generalize 
these  codes. 

Using  the  surface-hopping  model,  we  performed  a  comparative 
study  of  ion-pair  formation  on  two  similar  systems: 

Ar  +  I  Ar+  +  I„ 


(3) 


and 


K  +  I2  -  K  +  I2  (4) 

One  would  expect  these  two  reactions  to  be  rather  similar,  since  in 
each  case  the  single  active  electron  is  in  a  4s  orbital.  However,  an 
additional  final  channel  is  available  in  the  first  reaction.  Penning 
ionization  may  occur, 

*  + 

Ar  +  I„  -  Ar  +  I  +  e 
2  2 

Data  are  available  for  both  reactions.  The  FOM  group  in  Amsterdam  con¬ 
sidered  K  and  Dr.  K.  T.  Gillen  at  SRI  has  investigated  Ar  +  1^. 

In  a  collaborative  effort  with  Dr.  Gillen,  we  were  able  to  perform 
calculations  that  explained  some  of  the  differences  in  the  differential 
cross  sections  for  these  two  systems.  We  explored  how  differences  in 
the  matrix  elements  arising  from  the  inner  shell  vacancy  and  the  com¬ 
peting  ionization  channel  would  affect  the  cross  sections.  This  work 
has  been  submitted  to  the  Journal  of  Chemical  Physics,  and  a  preprint 
is  included  as  Appendix  C. 

(b)  Quantum  Mechanical  Perturbation  Theory  of  Nonadiabatic  Collisions 

Considerable  theoretical  effort^  has  been  directed  to  the 
quantum  mechanical  formulation  of  the  nonadiabatic  (i.e.,  multiple- 
potential-surface)  collision  problem.  The  most  general  theory,  and  the 
least  tractable,  is  a  full  coup led -channel  expansion  in  rotational, 
vibrational,  and  electronic  eigenstates.  We  sought  an  alternative 
formulation  that  preserved  the  essential  features  of  the  multiple- 
surface  crossings,  but  exploited  approximations  that  can  legitimately 

be  made  at  the  high  collision  energies  of  interest  (E  ~  50-100  eV) . 

16 

The  approach  chosen  was  the  impact  parameter  approximation,  in  which 
a  rectilinear  classical  trajectory  is  assumed  for  the  translational 
motion,  leading  to  a  time -dependent  Schroedinger 's  equation  for  the 
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remaining  degrees  of  freedom.  We  discovered  that  model  potentials  of 
the  same  type  used  in  previous  studies  of  collisions  of  alkali  atoms 
and  diatomic  molecules  led  to  a  particularly  simple  form  of  the  final 
equations.  These  equations  were  solved  using  the  Magnus  approximation.1' 

We  will  briefly  summarize  the  results  of  our  calculations. 

We  performed  calculations  using  model  potentials  chosen  to  model  the 
reaction 

Cs  +  02  -  C2+  +  02-  (6) 

Initially,  isotropic  potentials  were  assumed;  that  is,  molecular  rota¬ 
tion  was  neglected.  The  results  exhibited  strong  oscillations,  corre- 

18 

sponding  qualitatively  to  those  observed  in  recent  measurements. 

These  oscillations  have  been  related  to  the  vibrational  motion  of  the 
molecular  0^  or  0^  during  the  time  of  the  collision.  We  also  treated 
anisotropic  potentials  using  the  Infinite-Order-Sudden  approximation. 

We  obtained  rotationally  averaged  cross  sections  between  specific 
vibrational  states.  The  oscillations  observed  in  the  isotropic  case 
persisted,  although  they  were  less  prominent. 

These  rotationally  averaged  calculations  also  enabled  us  to  examine 
the  effect  of  different  initial  vibrational  levels  v^  on  the  ion  pair 
formation.  The  model  predicted  that  the  oscillations  will  become  less 
pronounced  as  v^  is  increased.  Although  no  measurements  are  available 
to  test  this  result,  our  calculations  suggest  a  direction  for  future 
experiments. 

The  theory  and  calculations  are  described  in  greater  detail  in  a 
paper  that  has  been  submitted  to  the  Journal  of  Chemical  Physics.  A 
preprint  is  included  as  Appendix  D. 
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(c)  Time -Independent  Quantum  Mechanical  Theory  of  Nonadiabatic 
Collisions 


The  time -independent  theory  of  nonadiabatic  collisions  that 
we  are  developing  is  a  generalization  of  the  usual  IOS  approximation 
and  is  complementary  to  the  impact  parameter  approach  described  above. 
The  impact  parameter  treatment  assumes  a  straight  line  trajectory,  but 
explicitly  considers  the  internal  state  structure,  whereas  methods  based 
on  the  sudden  approximation  properly  treat  the  dynamics,  but  assume 
internal  energies  are  degenerate. 


The  sudden  approximation  is  conceptually  appealing  for  our 

* 

application  to  systems  like  He  +  0^,  since  in  the  experiments  rotational 

states  of  the  target  are  not  resolved  and  in  the  sudden  approximation 

they  do  not  have  to  be  considered  individually.  Although  the  IOS 
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formalism  has  been  given  *  for  the  case  where  vibrational  channels 
are  considered  explicitly,  most  calculations  have  been  done  in  the  rigid 
rotor  approximation  (i.e.,  ignoring  vibration). 


In  the  present  work,  we  develop  the  time-independent  close 
coupled  equations  for  the  case  of  two  electronic  states  in  the  sudden 
approximation  for  all  internal  degrees  of  freedom,  that  is,  vibrational 
and  rotational.  We  start  with  a  Hamiltonian  expressed  in  terms  of  the 
matrix  elements  between  the  electronic  eigenstates.  We  then  seek  a 
solution  to  the  time-independent  Schroedinger  equation.  The  solution 


is  written  as  the  column  vector 


♦  (R.r) 
^(R.r)  » 


where  the  wavef unctions  on  the 


covalent  and  ionic  surfaces,  respectively,  are  given  by  expansions  in 


terms  of  the  internal  states  of  0^  and  0^ 


t0<Rir> 


L  F.(St)  cp.(r) 
i  1  i 


(7) 


_  Z 


7  G<  (R)x,<0 


MV) 


(8) 


In  this  case  c.  (r)  and  v  (r)  refer  to  all  internal  degrees 
i  j 

of  freedom;  that  is,  r  refers  to  all  internal  coordinates.  In  the 
sudden  approximation,  we  solve  the  simplified  set  of  two  coupled 
equations : 


f d2 

«■ 

1(1  +  1)  --2 

2u 

'VR’E)  voi(R’r)‘ 

’f(H.r)' 

dR2 

- —  +  k 

R^ 

m 

1  - 

.V01(R'r)  Vll(R»r). 

g(R,r) 

*  m 

0  (9) 


Here  i  is  an  appropriately  chosen  value  of  the  orbital  angular  momentum 
—  2  2P-  — 

and  k  »  — *r  (E-7) «  where  *  is  an  average  internal  energy  for  the  two 
A 

electronic  states,  for  the  functions  f(R,r)  and  g(R,r).  Replacing  & 
by  l  and  the  internal  energies  by  an  average  value  are  the  heart  of  the 
sudden  approximation.  In  Eq.  (9),  the  internal  coordinates  denoted  by 
r  are  considered  as  fixed  parameters.  The  desired  solutions  of  the 
original  problem  are  then  given  as: 


Fin(R)  *  <  ^  I  f  (R»r)  Cr )  > 


(10) 

Gjn(R)  -  <  Xj  (r>  1  S(R,r)  |cpn(r>  >  <u> 

Roughly  speaking,  Eqs.  (10)  and  (11)  average  the  solutions  f(R,r)  and 
g(R,r)  for  fixed  values  of  r  over  the  internal  wavefunctions  in  r.  In 
the  present  formulation,  the  number  of  coupled  equations  is  just  the 
number  of  electronic  surfaces  in  the  problem. 

The  S  matrix  elements  are  written  on  terms  of  four  N  x  N  blocks 
whose  elements  are  S„,  S^,  S_,  and  S^,  where  the  superscripts 
denote  the  initial  and  final  electronic  surface,  and  the  subscripts 
denote  the  initial  and  final  target  state.  In  terms  of  the  2  x  2  S 
matrix  calculated  by  solving  Eq.  (9),  whose  elements  are  Sqq,  S^, 

S^,  we  have 
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(12) 


<*i 

(r) 

■  soo 

1  Vj 

(r)) 

<9t 

(r) 

1S01 

'  Xj 

(r)> 

(r) 

'S10 

1  *J 

(r)> 

(r) 

lsll 

1  Xj 

(r)> 

In  this  work  we  are  interested  in  the  total  cross  section  for 
scattering  from  a  particular  vibrational  state  i  on  one  electronic  sur¬ 
face  to  vibrational  state  j  on  the  other  electronic  potential  surface, 
where  the  cross  section  has  been  summed  over  all  final  rotational  states. 
We-  now  note  that  the  internal  coordinate  stands  for  the  vibrational 
coordinate  as  well  as  the  body  fixed  angle,  y,  between  the  diatomic 
axis  and  the  vector  from  the  diatomic  center  of  mass  to  the  third 
body.  Then  the  desired  cross  section  is 


a1Q  (j  ~  i)  *  %  t?10(j  ~i;y)  d  cos  y 


where 


(j  -  i ; y)  -  — ff  E  (2£  +  1) 


e10  \  |2 
^ij  I 


Here  we  have  explicitly  written  the  dependence  of  S  on  the  angular 
momentum  i  and  angle  y. 

For  our  first  application  of  this  formalism,  we  used  the  following 
model  potential  functions: 

V00  (R>  -  CL  exp  (-C2R)  (18) 

V11  (R)  "  voo  (R)  "  1/R  +  AE  (19) 


v0i  00  -  C3  exp  ,-C4R)  sin  2y 


■13 


where  is  the  asymptotic  separation  between  the  electronic  states. 

The  parameters  *  319.8680  a.u.,  =  1.9319  a.u.,  =  0.35  a.u., 

C,  »  0.5  a.u.,  and  &E  =  0.12578  a.u.,  were  chosen  to  model  the  reaction 
4 

Cs  +  02  _  Cs+  +  o’  (21) 

The  calculations  now  completed  have  assumed  that  sin  2y  is  fixed  at  1.0. 
The  cross  sections  were  summed  explicitly  over  final  vibrational  states. 
Calculations  were  performed  for  a  system  with  the  reduced  mass  of 
Li  +  0^  to  reduce  the  amount  of  computation  required.  The  internal 
energy  f ^  was  taken  as  the  vibrational  energy  of  state  i  of  0^,  and 
the  internal  energy  of  c  was  taken  as  the  internal  energy  of  state  j  of 

°2”  * 

Results  are  shown  in  Fig.  1  for  initial  vibrational  states  0  and  1. 
The  energy  scale  on  the  top  of  the  figure  also  shows  the  total  energy 
for  the  Cs  +  0^  system  with  the  same  relative  velocity  as  that  of  the 
Li  +  mass.  Calculations  must  be  performed  at  higher  total  energies 
before  the  results  can  be  evaluated.  The  increase  of  the  cross  section 
with  initial  vibrational  state  is  also  apparent  in  the  figure.  Future 
work  also  includes  study  of  the  effect  of  the  choice  of  e  on  the  results 
and  the  evaluation  of  the  charge  exchange  cross  section  for  anisotropic 
potentials. 


14 


CHARGE  EXCHANGE  CROSS-SECTION  (a„ 


List  of  Papers  Completed.  1975-1980 


1.  F.  T.  Smith,  D.  L.  Huestis,  D.  Mukherjee,  and  W.  H.  Miller, 
"Semiclassical  Perturbation  Scattering  by  a  Rigid  Dipole," 
Phys.  Rev.  Lett.  35,  1073  (1975). 


2.  T.  F.  Gallagher,  R.  E.  Olson,  S.  A.  Edelstein,  and  R.  M.  Hill, 
"Collisional  Angular  Momentum  Mixing  of  Highly  Excited  Na  (n^D) 
by  N2  and  CO,"  Phys.  Rev.  A  16,  441  (1977). 


3. 


R.  E.  Olson,  "Theoretical  Excitation  Transfer  Cross  Sections 
for  Rydberg  Na  (n^D  n^F)  Transitions  from  Collisions  with  He, 
Ne,  Ar,"  Phys.  Rev.  A  15,  631  (1977). 


4.  A.  P.  Hickman  and  H.  Morgner,  "Calculation  of  Associative  and 
Penning  Ionization  of  H  and  D  by  He  (2^S)  and  He  (2^S)," 

J.  Chem.  Phys.  67,  5484  (1977). 


5.  R.  P.  Saxon,  R.  E.  Olson,  and  B.  Liu,  "Ab  Initio  Calculations 
for  the  X^E,  A^tt,  B^E  States  of  Na  Ar:  Emission  Spectra  and 
Cross  Sections  for  Fine -Structure  Transitions  in  Na-Ar  Collisions," 
J.  Chem.  Phys.  67,  2692  (1977). 

6.  W.  H.  Miller  and  F.  T.  Smith,  "Semiclassical  Perturbation  Theory 
of  Electron-Molecule  Collisions,"  Phys.  Rev.  A  JL7,  939  (1978). 

7.  D.  Mukherjee  and  F.  T.  Smith,  "Semiclassical  Perturbation  Scattering 
Theory  of  Electron-Polar  Molecule  Collisions:  Total  Excitation  and 
Scattering  Cross  Sections,"  Phys.  Rev.  A  JJ,  954  (1978). 

B,  A.  P.  Hickman  and  F.  T.  Smith,  '^omentum  Transfer  in  Electron- 
Polar  Molecule  Collisions:  Results  of  Semiclassical  Perturbation 
Scattering  Theory,"  Phys.  Rev.  A  17,  969  (1978). 

B.  Liu  and  R.  E.  Olson,  "Potential  Energies  for  Ca~+:  Cross 
Sections  for  Collisions  of  Cai"  and  Rydberg  Ca  with  Ca,  "  Phys. 

Rev.  A  18,  2498  (1978). 


R.  P.  Saxon,  K.  Kirby,  and  B.  Liu,  "Ab  Initio  Configuration  Inter¬ 
action  Study  of  The  Low-Lying  Electronic  States  of  MgH,"  J.  Chem. 
Phys,  69,  5301  (1978). 

A.  P.  Hickman,  "Theory  of  Angular  Momentum  Mixing  in  Rydberg  Atom- 
Rare  Gas  Collisions,"  Phys.  Rev.  A  18,  1339  (1978). 


12.  A.  P.  Hickman,  "Relation  Between  Low-energy  Electron  Scattering  and 
2-Changing  Collisions  of  Rydberg  Atoms,"  Phys.  Rev.  A  19,  994  (1979). 

13.  R.  E.  Olson,  "Rydberg  Atom-Rydberg  Atom  Ionization  Sections," 

J.  Phys.  B  12,  L109  (1979). 

14.  K.  Kirby,  R.  P.  Saxon,  and  3.  Liu,  "Oscillator  Strengths  and 
Photodissociation  Cross  Sections  in  the  y  —  A  and  y  —  B  Band  Systems 
in  MgH;  Astrophys.  J.  231,  637  (1979). 

15.  K.  Kirby,  W.  Roberge,  R.  P.  Saxon,  and  B.  Liu,  "Fhotodissociation 
Cross  Sections  and  Rates  for  CH"1"  in  Interstellar  Clouds,"  Astrophys.  J 
in  -press. 

16.  K.  Kirby,  R.  P.  Saxon,  and  B.  Liu,  "Excited  States  of  CH+:  Potential 
Curves  and  Transition  Moments,"  J.  Chem.  Phys.,  in  press. 

17.  A.  P.  Hickman,  "Approximate  Scaling  Formula  for  Angular  Momentum 
Changing  Collisions  of  Rydberg  Atoms,"  submitted  to  Phys.  Rev.  A. 

18.  A.  P.  Hickman  and  K.  T.  Gillen,  "Comparison  of  Ion-Pair  Formation 
in  the  Systems  Ar*  +  and  K  +  I2,"  submitted  to  J.  Chem.  Phys. 

19.  A.  P.  Hickman,  "Model  for  Fast,  Nonadiabatic  Collisions  Between 
Alkali  Atoms  and  Diatomic  Molecules,"  submitted  to  J.  Chem.  Phys. 


Ill  REFERENCES 


1.  T.  F.  Gallagher,  S.  A.  Edelstein,  and  R.  M.  Hill,  Fhys.  Rev.  Lett. 
35,  654  (1975);  Phys.  Rev.  A  15,  1945  (1977). 

2.  R.  E.  Olson,  Phys.  Rev.  A  _15,  631  (1977). 

3.  A.  P.  Hickman,  Phys.  Rev.  A  _18,  1339  (1978). 

4.  A.  P.  Hickman,  Phys.  Rev.  A  _19,  994  (1979). 

5.  X.  F.  Gallagher,  R.  E.  Olson,  W.  E.  Cooke,  S.  A.  Edelstein,  and 
R.  M.  Hill,  Phys.  Rev.  A  16,  441  (1977). 

6.  K.  A.  Safinya  and  T.  F.  Gallagher,  Phys.  Rev.  A.,  in  press. 

7.  K.  Kirby,  W.  Roberge,  R.  P.  Saxon,  and  B.  Liu,  Astrophys  J.,  in 

press . 

8.  F.  Biraben,  Xh4se,  Paris  VI  (1977). 

9.  A.  W.  Kleyn,  M.  M.  Hubers,  and  J.  Los,  Chem.  Phys.  34,  55  (1978). 

10.  J.  A.  Aten.  G.E.H.  Lanting,  and  J.  Los,  Chem.  Phys.  JL9,  241  (1977). 

11.  J.  A.  Aten  and  J.  A.  Los,  Chem.  Phys.  £5,  47  (1977). 

12.  G.  A.  Parker  and  R.  X.  Pack,  J.  Chem.  Phys.  68,  1585  (1978). 

13.  D.  Secrest,  J.  Chem.  Phys.  62,  710  (1975). 

14.  D.  J.  Kouri,  in  Atom-Molecule  Collision  Xheorv.  R.  B.  Bernstein,  ed. 
(Plenum,  NY,  1979). 

15.  For  example,  see  review  by  M.  S.  Child,  in  Atomic -Molecule  Collision 
Xheorv .  R.  B.  Bernstein,  ed.  (Plenum,  NY,  1979). 

16.  S.  Geltman,  logics  in  Atomic  Collision  Xheorv  (Academic  Press,  NY, 
1969). 

17.  P.  Pechukas  and  J.  C.  Light,  J.  Chem.  Phys.  44,  3897  (1966). 

18.  A.  W.  Kleyn,  M.  M.  Hubers,  and  J.  Los,  Chem.  Phys.  34,  55  (1978). 


18 


APPENDIX  A 


APPROXIMATE  SCALING  FORMULA  FOR  COLLISIONAL  ANGULAR 
MOMENTUM  MIXING  OF  RYDBERG  ATOMS 

A.  P.  Hickman 

Molecular  Physics  Laboratory 
SRI  International,  Menlo  Park,  CA  94025 

ABSTRACT 

An  approximate  scaling  formula  has  been  determined  that  permits 
the  rapid  estimation  of  cross  sections  for  angular  momentum  changing 
collisions  of  Rydberg  atoms  with  a  variety  of  targets,  using  informa 
tion  about  low  energy  electron  scattering  from  the  target.  The 
formula  is  obtained  by  fitting  the  results  of  coup led -channel  and 
3om  approximation  calculations  to  functions  of  reduced  parameters. 
Application  to  Na(nd)  +  He,  Ne,  Ar,  N^,  CH^,  C^H^;  Rb(nf)  +  He,  Ar, 
Xe;  and  Xe(nf)  +  CO^,  suggests  that  the  accuracy  is  about  a  factor 
of  two. 


MP  30-73 


I  INTRODUCTION 


Considerable  attention  has  been  devoted  recently  to  collisions 
involving  Rydberg  atoms,  especially  the  "i-mixing"  or  angular  momentum 

1-li 

changing  process 

R*  (ni)  +  X  -  R*  (n L*  )  +  X  L* it  (1) 

k 

where  R  is  the  Rydberg  atom  and  X  is  the  collision  partner.  It  has 
generally  been  argued  that  one  should  be  able  to  analyze  this  process 
in  terms  of  the  cross  section  for  low  energy  electron  scattering  from 
X.  In  this  paper  we  present  a  very  simple  scaling  formula  that  can  be 
used  to  estimate  the  cross  section  for  reaction  (1)  to  about  a  factor 
of  two,  starting  from  information  about  e-X  scattering.  This  degree 
of  accuracy  is  useful  because  the  cross  sections  may  vary  by  as  much 
as  an  order  of  magnitude  over  the  range  of  n  considered. 

The  scaling  formula  was  determined  by  fitting  the  results  of 
coupled -channel  and  Born  approximation  calculations  to  functions  of 
reduced  parameters.  These  parameters  have  a  reasonable  physical 
interpretation  in  terms  of  the  "nearly  free  electron"  picture.  The 
determination  of  the  formula  and  its  physical  interpretation  are  dis¬ 
cussed  in  Section  II.  Section  III  contains  the  results  of  calculations 
for  several  systems,  including  molecular  targets.  The  accuracy  of  the 
results  for  more  complicated  targets  suggests  that  their  Internal 
structure  may  play  only  a  minor  role  in  the  collisions  studied. 


II  THE  SCALING  FORMULA 


A.  Determination 

6,8 

Our  previous  calculations  of  angular  momentum  mixing  may  be 
characterized  as  exact  or  approximate  solutions  of  a  specific,  well- 
defined  model  problem.  We  will  begin  by  summarizing  this  problem,  and 
then  show  how  one  might  expect  to  estimate  the  desired  cross  sections 
using  empirical  functions  of  reduced  parameters.  Then  one  additional 
approximation  is  introduced  to  tie  the  model  problem  to  the  physical 
system  involving  the  Rydberg  atom  and  an  arbitrary  collision  partner. 
The  model  problem  may  be  posed  as  follows.  Consider  the  system 

it  it 

Na  +  X  indicated  in  Figure  1.  The  Rydberg  atom  Na  has  principle 
quantum  number  n,  and  we  consider  only  the  initial  level  2-2  and  other 
levels  2*3,4,...,n-l.  Assume 

Eni  ‘  End  “  *E’  >3,4,...,n-l  (2) 

This  is  approximately  true  for  sodium  since  the  d  quantum  defect  5.  is 

d 

much  larger  chan  5^  for  2  >  2.  The  interaction  potential  is  given  by  a 
Fermi  pseudopotential 

V(R,r )  -  2rtA  5(r-R)  (3) 

For  the  present  discussion,  A  is  to  be  considered  an  arbitrary  constant 
We  wish  to  calculate  .  ,  defined  by 

n-1 

•j.„u  ■  2  »<nd-'u>-  (4> 

2-3 
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using  che  quanttan  mechanical  Arthurs  and  Dalgarno  formalism  or  the 


Bom  approximation  as  discussed  previously 


6,3 


We  expect  that  a 


2-mix 


will  depend  on  a.  A,  aE,  and  v,  the  relative  velocity  of  Na  and  X. 
That  is, 

J,  -  a,  .  (n,A,  AE,v)  (3) 

A -mix  2 -mix 

We  wish  to  find  a  functional  form  for  o^_.  ♦  Previously,  we 
0 

showed  that  if  4E  ■  0,  the  results  of  the  Bom  approximation  may  be 
well-fitted  by 

/  \o 

.733 


a .  .  a 

2-mix 


(■] 


(aE-0) 


(6) 


We  can  write  this  in  a  slightly  different  form  by  separating  out  a 
4  2 

factor  rm  a^  (the  so-called  geometrical  cross  section  of  the  Rydberg 
atom) ,  and  making  A/v  dimensionless  by  dividing  by  che  unit  length  and 
velocity  in  atomic  units: 

<7. 


4  2  2 

a  TO  a„  3 
2-mxx  0 


(AE-0) 


(7) 


where 


m*  a'T  2  6.733 
e  0  v  n 


(8) 


The  constant  mg  is  che  electron  mass  in  atomic  units.  If  aE  t  0  our  pre- 

g 

vious  results  suggested  that  in  che  weak  coupling  limit  a  more  general 


formula  could  be  written 


where 


2  „ 
n  aQAE 

hv 


(10) 


f(Y)  is  an  approximation  to  the  family  of  functions  shown  in  Figure  2 
of  reference  8  that  depend  weakly  on  the  parameter  n.  In  neglecting 
the  dependence  on  n  we  simplify  the  functional  form  but  introduce  some 
ambiguity  into  the  determination  of  f.  Eq.  (9)  is  expected  to  be  a 
reasonable  approximation  in  the  weak-coupling  (small  A  or  large  n) 
limit.  We  now  postulate  that  a  more  general  formula  can  be  obtained 
by  writing 


Vmix  *  *^*1  S(3)f(Y)  (U) 

where 

S(3)  -  ki1  as  8-  0  (12) 

We  will  find  later  the  constant  y  has  the  value  0.715.  We  performed  a 
number  of  calculations  using  the  coupled  channel  method  and  Bom 
approximation  to  test  the  usefulness  of  Eq.  (11).  For  the  situations 
previously  considered,  the  Bom  approximation  could  be  applied  only 
when  n  was  large,  and  the  coupled  channel  method  was  feasible  only  when 
n  -and  the  number  of  channels  -was  small.  However,  by  varying  A  and  v, 
we  have  been  able  to  probe  a  wider  region  of  3-y  space  and  still  keep 
the  number  of  channels  small.  The  results  of  these  calculations  are 
listed  in  Table  1. 
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Is  was  possible  co  find,  empirically,  functions  f(y)  and  sO)  such 
chat  Che  results  of  the  coupled  channel  and  3om  approximation  calcula¬ 
tions  were  approximated  by  Eq.  (11)  to  about  a  factor  of  two.  These 
functions  are  shown  in  Figures  2  and  3  and  tabulated  in  Tables  2  and  3. 

At  this  point  we  have  obcained  an  approximate  scaling  rule  that 
can  be  used  co  estimate  the  results  of  a  model  problem  involving  a 
Fermi  pseudopocential  with  arbitrary  constant  A.  We  now  make  an  addi¬ 
tional  approximation  to  relate  A  to  the  low-energy  electron  scattering 
by  X.  Let 


i  ■> 

4  tt  A”  ■  3  ,  (l/2n“) 
el 


(13) 


a  .  is  the  e-X  scattering  cross  section  at  the  electron  energy  l/2n  , 
e  1 

which  is  the  average  kinetic  energy  of  the  electron  in  the  quantum  level 
n.  As  a  -  A  -  L,  the  scattering  length.  Some  justification  is  given 
for  Eq.  (13)  in  reference  8.  In  this  work  we  find,  a  posteriori,  that 
it  is  reasonably  successful. 

With  the  substitution  defined  by  Eq.  (13),  Che  final  form  of  the 
scaling  rule  is  given  by  Eqs.  (10),  (11),  and 


1 


m  a*  3.337 
e  0  vn 


4n 


(14) 


Although  the  scaling  rule  is  empirical,  it  can  be  related  to  a 
natural  physical  interpretation  of  the  scattering  process.  This 


interpretation  will  be  discussed  in  the  next  section.  Finally,  note 


that  although  the  preceeding  discussion  has  assumed  that  the  Rydberg 
atom  was  sodium,  we  find  that  the  formula  may  be  applied  to  other  atoms 
as  well.  The  essential  feature  of  the  model  is  that  the  initial  level 
is  separated  by  &E  from  a  nearly  degenerate  set  of  possible  final  levels. 
For  large  n,  the  difference  in  the  number  of  final  states  for  an  initial 
£«2  or  3  is  small.  Hence  we  will  test  the  formula  for  collisions  in¬ 
volving  Na  (nd),  Rb  (nf),  and  Xe  (nf). 

3.  Physical  Interpretation  and  Discussion 

The  approximate  scaling  formula  is  a  product  of  three  factors. 

4  2 

The  geometrical  factor  rrn  aQ  shows  that  the  cross  section  scales  with 
the  overall  sire  of  the  atom.  We  interpret  the  parameter  3  as  a 
coupling  strength,  and  g(3)  as  the  probability  that  the  collision 
partner  will  encounter  the  orbiting  Rydberg  electron.  3  increases 
with  which  gives  an  effective  size  of  the  collision  partner,  and 

decreases  with  n,  because  for  larger  n  the  electron  "cloud"  is  more 
diffuse.  We  interpret  the  function  f(y)  as  the  probability  that  an 
elastic  collision  between  the  orbiting  electron  and  the  collision 
partner  will  cause  a  transition  into  a  new  energy  level.  It  is  inter¬ 
esting  to  note  that  the  parameter  y  may  be  rewritten 

v 

V  ?  (A6>  7“  (15) 

where  <15  Is  the  difference  in  the  quantum  defects  of  initial  and  final 
states,  and  v  *  1/n  (atomic  units),  the  velocity  corresponding  to  the 

d 


average  kinetic  energy  of  an  electron  with  principle  quantum  number  n 


This  formula  is  obtained  by  expanding 

1  1 

AE  “  2(n+5^)2  '  2(n+5^  (16 

in  the  limit  «  n  and  substituting  into  Eq.  (10).  When  v  is  large, 

it 

f(y)  -0  and  hence  the  2-mixing  will  be  small.  Eq .  (15)  shows  that 
this  can  occur  because  the  inelasticity  is  large,  or  because  v^  is 
large  compared  to  v,  and  hence  it  is  less  likely  that  a  collision  will 
deflect  the  electron  enough  to  change  the  shape  of  its  orbit.  Con¬ 
versely,  the  2-mixing  will  be  larger  when  y  is  small.  This  may  occur 
either  because  the  energy  difference  be ween  initial  and  final  states 
is  small,  or  because  the  electron  is  moving  3lowly  relative  to  the 
collision  partner  and  a  collision  tends  seriously  to  perturb  its 
orbit. 

It  is  instructive  to  consider  various  limiting  values  of  the 

scaling  functions  g  and  f.  We  have  already  noted  that  the  function  g 

is  parabolic  as  the  argument  approaches  zero.  For  large  values  of 

the  argument,  g  tends  to  saturate  at  a  value  of  0.5  Co  0.6.  This 

4  2 

number  is  somewhat  arbitrary  since  the  prefactor  rrn  a^  could  equally 
4  2  4  2 

well  have  been  2rm  aQ  or  4ttu  a^.  Appropriate  limiting  values  of  f(y) 
are  more  easily  defined.  As  iE  -  0  (or  v  -  *) ,  y  -  0,  f(y)  —  1.  We 
can  thus  draw  the  following  conclusions  about  the  general  behavior  of 
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1-mixing  cross  sections 4  For  large  n,  £  —  1,  and  the  angular  momentum 
levels  nd  and  ni'  (4*  a  2_)  are  defectively  degenerate.  The  decrease 
of  <3  .  (n)  with  increasing  n  is  due  to  the  reduced  coupling  strength 

of  the  diffuse  electron  cloud.  On  the  other  hand,  the  (inelastic) 
cross  section  tends  to  be  small  for  small  n  because  of  the  smaller 
geometrical  cross  section  and  the  increased  importance  of  the  energy 
level  splitting  &E. 

We  now  consider  how  the  velocity  dependence  of  the  4-mixing  cross 

section  is  controlled  by  the  functions  £  and  g.  At  large  n,  y  will 

generally  be  sufficiently  small  so  chat  f  a»  1  for  typical  (thermal) 

-3 

value  of  v.  (Note  that  «  n  ).  Then  the  behavior  of  g(3)  shows 

that  decreases  for  larger  values  of  v.  In  the  limit  v  - 

2 

a,-mix  *  ^'V  buC  ^°r  sma^er  va^ue*  of  v  the  dependence  may  go  as  1/v 
or  weaker.  The  1/v  dependence  corresponds  to  the  intuitive  notion 
that  a  slower  projectile  spends  more  time  passing  through  the  Rydberg 
atom,  and  consequently  has  a  higher  probability  of  encountering  the 
electron.  The  opposite  may  be  true  for  small  n  when  the  coupling  is 
stronger.  In  this  case,  g(@)  may  achieve  is  saturated  value  of 
—  0.5-0. 6  for  a  range  of  thermal  velocities,  so  chat  the  velocity 
dependence  of  gj  ^  will  be  determined  by  f (y) .  Examination  of 
Figure  (  )  shows  that  will  then  increase  with  larger  values  of 

v,  because  increasing  v  has  the  same  effect  as  decreasing  &E.  Similar 
behavior  has  been  analyzed  theoretically  in  other  inelastic  collisions 


by  Olson 
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Ill  RESULTS  AND  DISCUSSION 


from  He  and  there  is  not  a  strong  energy  dependence  of  the  cross 

section.  For  e-CO^  and  e-CH, ,  we  used  the  scattering  data  compiled  by 
16 

Itikawa.  These  data  are  the  total  momentum  transfer  cross  section, 
which  may  include  Inelastic  processes,  rather  than  the  elastic  cross 
sections  required  by  the  theory.  For  the  CH^  data,  Itikawa  estimates 
an  uncertainty  in  the  data  of  about  a  factor  of  two.  At  this  level  of 
accuracy,  and  at  the  very  low  electron  energies  involved,  the  substitu¬ 
tion  of  the  momentum  transfer  cross  section  for  the  elastic  cross 
section  is  probably  not  too  serious  an  approximation.  This  view  is 
supported  by  Itikawa.  The  CO^  data  is  more  accurate,  but  the  un¬ 
certainty  in  the  scaling  formula  probably  renders  the  distinction 
between  momentum  transfer  cross  section  and  elastic  cross  section 
irrelevant.  Finally,  the  e  -  C  H  scattering  data  of  McCorkle  et  al.^ 

J  0 

was  used.  In  this  case  also,  the  momentum  transfer  cross  sections  were 
measured . 

The  energy  differences  &E  of  Ha  were  taken  to  be  the  d-f  splittings 

18 

measured  by  Gallagher  et  al.  The  parameter  y  in  the  Rb  and  Xe  colli¬ 
sions  were  determined  from  Eq .  (15),  assuming 

Afi  *8  5f  (18) 

where  5^  ■  0.02  for  Rb  and  0.055  for  Xe. 
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3.  Cross  Sections 


Results  for  collisions  of  Na  and  Rb  Rydberg  atoms  with  rare  gases 
are  presented  in  Figures  4-9.  The  agreement  for  the  case  that  the 
collision  partner  is  He  is  especially  good.  Xote  chat  the  present 
results  for  Ha(nd)  +  He  are  determined  using  the  scaling  formula; 
they  differ  only  slightly  from  the  coupled  channel  and  Born  approxi¬ 
mation  results  presented  previously.  The  good  agreement  for  He 
suggests  chat  Che  Fermi  pseudo potential  is  a  rather  good  approximation 
when  the  low  energy  electron  scattering  has  a  weak  dependence  on 
energy.  Also,  the  small  polarizability  of  He  supports  the  use  of  a 
short  range  (delta  function)  approximation  to  the  e-He  potential. 

For  Che  collisions  involving  Ar  and  Xe  the  energy  dependence  of  the 
electron  scattering  cross  section  does  Influence  the  predicted  result 
for  n  —  10-20.  The  calculated  cross  sections  are  in  reasonable  agree¬ 
ment  with  Che  data  in  this  range,  suggesting  chat  the  approximate 
method  of  including  the  energy  dependence  [Eq.  (13)1  is  at  least 
qualitatively  correct. 

We  have  also  applied  the  scalit g  formula  to  collisions  involving 
more  complicated  partners.  In  these  calculations  Che  .ntemal  structure 
of  the  target  is  neglected,  although  it  may  of  course  influence  the 
electron  scattering  cross  sections  used.  The  qualitative  agreement 


observed  between  the  calculations  and  experiment  tends  to  suggest 


that  the  elastic  scattering  of  the  electron  by  the  target  is  the 
dominant  mechanism  of  2-mixing  in  these  collisions  as  well  as  in  chose 
involving  rare  gases.  Note  that  the  formula  is  not  intended  to  be 
used  for  a  collision  partner  that  has  a  charge  or  large  dipole  moment. 

Figures  10  thru  13  show  the  calculated  results  for  a  number  of 
molecular  targets.  The  agreement  with  experiment  is  generally 
reasonable.  The  excellent  results  for  Na(nd)  +  tend  to  confirm, 
as  in  the  case  of  He,  the  reliability  of  the  Fermi  pseudopotential 
when  the  low  energy  electron  scattering  does  not  have  a  strong  energy 
dependence . 

It  is  interesting  to  consider  the  structure  observed  in  the  experi¬ 
mental  dace  near  n  »  15  for  tfa(nd)  +  C,H  .  We  have  considered  what 

3  o 

form  of  the  elastic  cross  section  <y  would  be  necessary  to  lead  to 

the  observed  form  of  a,  .  (n),  assuming  the  validity  of  Eq.  (11). 

We  found  chat  the  unusual  structure  in  a.  .  (n)  for  n  —  13-16  could 

2 -mix 

be  fit  by  assuming  an  elec cron -pro pane  elastic  scattering  cross  section 

that  exhibits  a  strong  change  of  slope,  but  not  a  dip,  at  an  electron 

energy  0.060  eV,  which  is  the  average  electron  kinetic  energy  for 

n  •  15.  This  effective  elastic  cross  section  is  about  a  factor  of 

two  less  than  u  of  McCorkle  ec  al.^7  for  E  >  0.060  eV,  but  rises 
tn 

more  sharply  for  E  <  0.060  e?.  In  this  region  (n  —  15-17)  the  size 
of  the  coupling  strength  parameter  3  is  sensitive  to  two  competing 


factors.  As  n  increases,  A  increases  because  of  che  rapid  increase 
in  j  as  the  electron  energy  goes  to  zero.  However,  this  is  nearly 

ft  1 

-3 .37 

counteracted  by  the  factor  n  ’  that  reflects  the  weaker  effect  of 
the  increasingly  more  diffuse  electron  cloud.  The  net  effect  is  that 
the  coupling  strength  teeters  between  these  competing  influences. 

The  results  of  this  analysis  are  only  qualitative,  of  course,  but  they 
indicate  that  the  unusual  structure  in  the  Rydberg  cross  section  can  be 
related  co  a  plausible  behavior  of  the  corresponding  electron  scattering 
cross  section. 
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IV  CON  CUJ  DING  REMARKS 


An  approximate  scaling  formula  for  collislonal  angular  momentum 
mixing  of  low-4  Rydberg  atoms  has  been  determined  that  gives  reason¬ 
able  results  for  a  wide  variety  of  systems.  Although  the  formula  was 
determined  using  calculations  that  assumed  the  collision  partner  was 
a  rare  gas,  qualitative  agreement  is  also  obtained  for  targets  with 
internal  structure.  All  the  cross  sections  exhibit  a  maximum  as  a 
function  of  n,  but  the  shape  and  position  of  the  peak  may  vary  con¬ 
siderably.  For  example,  Rb  +  He  has  a  sharp  peak  at  n  *  11,  while 
Xe  +  CO ^  has  a  very  broad  maximum  for  cv»28.  The  formulas  obtained 
here  show  how  such  contrasting  behavior  may  be  qualitatively  related 
to  the  energy  level  splittings  of  the  Rydberg  atom,  the  energy 
dependence  of  the  electron  scattering  from  the  collision  partner, 
and  the  relative  velocity  of  the  collision. 
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n 

4 

5 

6 

7 

8 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 


6.867x10 
6.867x10 
6.867x10 
6.867x10 
6.867x10 
6.867x10 
6.867x10 
6.867x10 
6.867x10 
1.373x10 
4.856x10 
3.140x10 
3.140x10 
3.140x10 
3.140x10 
3.140x10 
4.856x10 
1. 717x10 ' 
1.030x10 


-4 

-4 

-4 

-4 

-4 

-4 

-4 

-4 

-4 

-3 

-4 

-4 

-4 

-4 

-4 

-4 

-4 

-3 

-3 


Lfa.u.'! 


1.19 
1.19 
1.19 
1.19 
1.19 
0.595 
0.3 
0.15 
0.075 
1.19 
0.841 
0.544 
0.272 
0.137 
0 .0686 
0.0343 
0.421 
1.488 
0.892 


I 

l 


TA3LE  1 


iction  of  n,  v. 


and  L  obtained  by  coup led - 


V 

3 

.  /t» 

— i-aix  — 

4.75 

16.21 

0.192 

3.82 

7.64 

0.387 

3.20 

4.14 

0.440 

2.75 

2.46 

0.402 

2.41 

1.56 

0.352 

3.20 

2.07 

0.324 

3.20 

1.05 

0.200 

3.20 

0.523 

0.082 

3.20 

0.260 

0.025 

1.60 

2.07 

0.535 

4.53 

4.14 

0.217 

7.00 

4.14 

0.102 

7.00 

2.07 

0.085 

7.00 

1.04 

0.055 

7.00 

0.523 

0.024 

7.00 

0.260 

0.0078 

4.53 

2.07 

0.171 

1.28 

2.07 

0.567 

2.14 

2.07 

0.462 
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TABLE  2 


Values  of  Che  collision  screngch  g  as  a  function  of  3*  For 
2 

3  £  0.5,  3  =  0.715  3”.  For  3  a  3.0,  we  arbitrarily  sec  g  ■  0.60 


a 

jh 

£ 

0 

0 

0.5 

0.18 

1.0 

0.33 

1.5 

0.43 

2.0 

0.51 

2.5 

0.56 

3.0 

0.60 
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TABLE  3 


Values  of  the  collision  efficiency. 


n  a  &E/ftv 

f 

0 

1.00 

1 

0.99 

2 

0.91 

3 

0.78 

4 

0.65 

3 

0.54 

6 

0.45 

7 

0.38 

3 

0.32 

9 

0.28 

10 

0.25 

11 

0.22 

12 

0.19 

13 

0.18 

14 

0.17 
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FIGURE  CAPTIONS 


1.  Schematic  diagram  of  coordinates  used  to  describe  the  Rydberg 
atom  and  its  collision  partner. 

2.  The  function  f(v),  which  we  interpret  as  the  probability  that 

a  collision  between  the  electron  and  the  collision  partner  will 
cause  a  change  of  angular  momentum  level  of  the  electron. 

3.  The  function  g(3)»  which  we  interpret  as  the  probability  of 
encounter  between  the  Rydberg  electron  and  a  collision  partner. 

4.  The  4-mixing  cross  sections  for  Na(nd)  +  He.  The  calculation 

of  de  Prune l4  and  Pascale  obtained  upper  and  lower  limits  to  the 
cross  sections. 

5.  The  4 -mixing  cross  sections  for  Na(nd)  +  Ne.  The  calculation 
of  de  Prune le  and  Pascaie  obtained  upper  and  lower  limits  to 
the  cross  sections. 

6.  The  4-mixing  cross  sections  for  Na(nd)  +  Ar.  The  calculation 

of  de  Prune l4  and  Pascale  obtained  upper  and  lower  limits  to  the 
cross  section. 

7.  The  4-mixing  cross  sections  for  Rb(nf)  +  He.  The  calculation 
of  de  Prune l4  and  Pascale  obtained  upper  and  lower  limits  to 
che  cross  sections. 

8.  The  4-mixing  cross  sections  for  Rb(nf)  +  Ar.  The  calculation 
of  de  Prune 14  and  Pascale  obtained  upper  and  lower  limits  to 

i  Che  cross  sections. 

9.  The  4-mixing  cross  sections  for  Rb(nf)  +  Xe.  The  calculation 

'  of  de  Prunel4  and  Pascale  obtained  upper  and  lower  limits  to 

the  cross  sections. 

10.  The  4-mixing  cross  sections  for  Na(nd)  +  N2 .  The  calculation 
of  de  Prune l4  and  Pascale  obtained  upper  and  lower  limits  to 
the  cross  sections. 

11.  The  4-mixing  cross  sections  for  Na(nd)  +  CH^. 
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12.  The  3-mixing  cross  sections  for  Na(nd)  +  C^H^.  Tiie  anusual  dip 
at  n*15  is  discussed  in  the  text. 

13.  The  3-mixing  cross  sections  for  Xe(nf)  +  CO,,. 


I 
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‘  i 
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Rb  (nf)  +  Xe 
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Abstract 


Wavefunctions  and  potential  curves  for  the  ground  X  I  state  and 
eight  excited  states  of  and  *"  symmetry  of  CH+  have  been  obtained 
using  ab  initio  conf iguration  interaction  (Cl)  methods.  In  order  to 
take  proper  account  of  valence-Rydberg  mixing,  Rydberg  functions  were 
included  in  the  basis  set.  The  orbitals  used  in  the  set  of  reference 
configurations  for  the  Cl  wavefunctions  included  both  valence  and  Rydberg 
orbitals,  determined  from  a  multiconfiguration  self-consistent  field  and 
a  natural  orbital  calculation,  respectively.  Transition  moments  between 
the  ground  electronic  state  and  the  2*E+,  3*E+,  A*!  and  2*11  states  have 
been  computed,  and  the  importance  of  these  states  in  photodissociation 
of  CH  is  discussed. 
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I .  Introduction 

The  CH+  molecule  has  been  of  great  interest  to  astrophysicists 

ever  since  its  discovery  in  diffuse  interstellar  clouds.*  In  particular, 

the  abundance  of  CH  appears  to  be  as  much  as  30  times  larger  than 

predictions  based  on  probable  formation  and  destruction  mechanisms  in 

2 

interstellar  clouds.  Also,  many  absorption  lines  observed  in  diffuse 

3 

interstellar  clouds  remain  unidentified.  Transitions  to  excited  states 

of  CH+,  not  yet  discovered  in  the  laboratory,  have  been  suggested  as 

4 

possible  origins  of  some  of  these  lines.  A  detailed  understanding  of 
the  structure  and  radiative  transitions  of  this  molecule  could  lead  to 
insights  regarding  these  astrophysical  observations. 

The  ground  electronic  state  of  CH+  has  *£+  symmetry  and  the  oscill¬ 
ator  strength  of  the  only  known  dipole  transition  A*H  -  X*£+  has  been 
5  6 

calculated  and  measured.  Accurate  potential  curves  for  the  four  lowest- 

lying  electronic  states,  X*I+,  A*H,  ^2  and  Z+,  dissociating  to  the  lowest 

+  22  7 

separated-atom  limit  C  (  P)  +  H(  S)  have  been  calculated  by  Green  et  al. 

No  higher-lying  excited  states  are  reliably  known,  although  Lorquet  et  al.^ 

produced  some  excited  state  potential  curves  showing  a  bound  *1  state 

1  +  9 

dissociating  to  C(  D)  +  H  .  Watson,  Stewart  and  Dalgamo,  using  the 

random  phase  approximation,  obtained  large  oscillator  strengths  for 

transitions  in  CH+  from  the  X*I+  state  to  higher  excited  states.  However, 


the  calculation  was  carried  out  at  only  one  intemuclear  separation,  so 


•'  T 


that  che  repulsive  or  attractive  nature  of  these  excited  states  was  not 
determined. 

In  this  work,  we  have  undertaken  a  study  to  characterize  those 
excited  states  of  CH+,  which,  in  dipole  transitions  from  the  ground  state, 
can  give  rise  to  band  spectra  or  to  direct  dissociation  of  the  molecule. 

We  have  calculated  potential  curves,  wavefunctions  and  transition  moments 
for  all  eight  states  of  ^£+  or  ^2  symmetry  lying  within  ~17  eV  of  the 
ground  state.  In  Section  II  we  discuss  the  basis  set  choice  and  the 
construction  of  configuration  interaction  wavefunctions  which  allow  for 
adequate  description  of  valence  and  Rydberg  state  interaction.  Additional, 
larger  calculations  designed  to  test  the  results  obtaindd  are  also  des¬ 
cribed.  The  results  of  all  the  calculations  as  well  as  the  astrophysical 
implications  are  detailed  in  Section  III. 

II .  Calculations 

The  CH+  excited  states  of  interest  in  studies  of  the  interstellar 

medium  must  lie  within  13.6  eV  of  the  ground  state  because  there  are  no 

photons  in  the  interstellar  medium  with  hv  >  13.6  eV  due  to  che  ionization 

of  atomic  hydrogen.  However,  in  order  to  compute  such  states  accurately 

the  interaction  with  states  lying  somewhat  higher  in  energy  has  to  be 

considered.  Figure  1  shows  the  lowest  six  CH+  separated  atom  limits 

which  give  rise  to  states  of  singlet  £  or  H  symmetry  and  the  measured 

energy  differences.  The  lowest  two  Rydberg  limits  C(^P,  2p3s)  +  H  and 
+  2 

C  (  P)  +H(23,2p),  lie  within  0.2  eV  of  each  other  and  within  about 


10  eV  of  the  ground  state  asymptote.  States  arising  from  these  Rydberg 
asymptotes  may  be  of  astrophysical  interest.  In  addition,  states  arising 
from  the  lower-lying  asymptotes  probably  experience  some  valence -Rydberg 
mixing -at  small  intemuclear  separations.  Thus,  it  is  necessary  to  design 
a  calculation  which  gives  a  balanced  treatment  of  both  valence  and  Rydberg 
states. 

Electronic  energies  and  wavefunctions  were  calculated  using  the 

configuration  interaction  (Cl)  method.  Each  wavefunction  was  expanded 

in  an  orthonormal,  n-particle  basis  set  of  C  symmetry-and  equivalence- 

csv 

restricted  configuration  state  functions  (CSF).  These  CSF ' s  were  linear 
combinations  of  Slater  determinants  which  have  the  appropriate  molecular 
symmetry  and  multiplicity.  The  Slater  determinants  were  constructed 
from  an  orthonormal  set  of  molecular  orbitals  which  were  expanded  in 
terms  of  an  elementary  basis  set  of  Slater-type  functions  (STF)  centered 
at  the  atonic  nuclei. 

The  extended  SCF  basis  set  used  in  these  calculations  is  given  in 

Table  1.  The  (5s/4p)  basis  of  Clementi  and  Roetti^  for  C(^P)  was 

augmented  by  two  3d  functions,  with  exponents  chosen  to  account  for  the 

distortion  due  to  both  long  range  electrostatic  interaction  and  short 

range  chemical  interactions.  Rydberg  3s  and  3p  functions  with  exponents 

1  3 

optimized  for  the  SCF  energies  of  carbon  atom  Rydberg  states  P  and  0, 
respectively,  were  added  to  the  set.  The  hydrogen  aton  basis  was  chosen 
to  describe  both  the  separated-atom  and  molecular-bonding  regions,  as 
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well  as  the  Rydberg  2s  and  2p  limits.  The  molecular  basis  set  consisted 
of  24a,  12rr  and  3§  functions. 

In  the  construction  of  the  Cl  wavefunctions ,  care  was  taken  to  treat 

Rydberg  and  valence  states  in  an  equivalent  manner  in  order  to  properly 

describe  valence-Rydberg  mixing.  The  configuration  lists  for  the  and 

symmetries  consisted  of  all  single  and  double  excitations  with  respect 

to  the  lists  of  reference  configurations  given  in  Table  2,  in  which  both 

valence  and  Rydberg  orbitals  were  occupied.  In  order  to  dissociate  the 

molecule  correctly,  in  these  configurations,  the  3 a  and  4 a  orbitals, 

which  asvmototically  have  the  form  3a  ■  2p  +  Is  and  4a  *  2p  -  Is  , 

c  H  c  H 

were  treated  equivalently.  The  5a,  6a,  2::  and  3rr  orbitals  are  Rydberg 

orbitals.  The  number  of  configurations  in  the  Cl  wavefunction  was  6808 

for  the  2  symmetry  and  9573  for  the  3  symmetry. 

For  each  intemuclear  separation,  the  orbitals  for  these  calculations 

were  determined  by  a  four-step  procedure  designed  to  provide  physically 

realistic  Rydberg  as  well  as  valence  orbitals.  In  the  first  step,  valence 

a  orbitals  were  determined  by  a  multiconfiguration  self-consistent  field 

(MCSCF)  calculation  on  the  X^2+  state  which  included  the  configurations 
2  2  2  2  2  2 

la  2a  3a  end  la  2a  4a  ,  required  for  proper  dissociation,  along  with  the 

2  2  2  2  2  2 

correlating  configurations  la  3a  Itt  and  la  4a  Itt  .  With  the  a  orbitals 

frozen  in  the  second  step  the  valence  n  orbital  was  determined  from  a 

properly  dissociating  MCSCF  calculation  on  the  A *11  state  which  included  the 
2  2  2  2 

configurations  la  2a  3alrr  and  la  2a  4a Itt.  In  the  third  step,  two  Rydberg 
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a  orbitals  were  obtained  as  natural  orbitals  from  a  frozen  core  two 


electron  Cl  calculation  for  4  Rydberg  states  which  included  all  config¬ 
urations  with  one  electron  in  the  3c  or  4c  orbital  and  one  electron  in  an 
external  c  orbital,  i.e.  no  for  n  2  5.  In  this  calculation,  six  electronic 
states  were  obtained  and  the  occupied  Rydberg  orbital  from  the  lowest  two 
states  in  which  the  Rydberg  orbital  was  occupied  was  selected;  i.e.,  one 
orbital  was  taken  from  each  state.  Finally,  two  Rydberg  tt  orbitals  were 
obtained  in  an  analogous  manner  from  a  Cl  calculation  for  Rydberg 
states  which  included  all  configurations  with  one  electron  in  the  3c  or 
4<j  orbital  and  one  electron  in  an  external  rr  orbital,  i.e.,  nrr  for  n  2  2. 

The  orbitals  thus  obtained  were  used  for  the  Cl  calculations  on  both 
the  2  and  H  symmetries.  The  calculations  were  performed  at  25  internuclear 

separations  between  R  *  1.5  a  and  20  a  . 

o  0 

A  limited  number  of  larger  configuration  interaction  calculations 

were  carried  out  to  test  whether  certain  improvements  which  could  be  made 

in  these  calculations  would  significantly  affect  the  results.  As  will  be 

discussed  in  the  next  section,  the  basis  set  used  resulted  in  errors  in 

asymptotic  excitation  energies  as  large  as  0.4  eV.  Therefore,  in  these 

extended  calculations  an  enlarged  basis  set  listed  in  Table  1  was  employed. 

10  3 

This  set  consisted  of  the  (6s/4p)  basis  of  Clementi  and  Roetti  for  C(  P), 

augmented  by  two  3d  functions,  two  3s  and  two  3p  Rydberg  functions  and 

one  4f  function.  The  diffuse  3d  function  was  chosen  to  maximize  the 

3 

polarizability  of  the  C(  P)  atom  and  the  tight  3d  function  to  maximize 
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the  energy  lowering  in  a  Cl  calculation  on  the  CH  molecule  (  Z  state). 


The  tighter  3s  and  3p  Rydberg  functions  were  optomized  with  respect  to 
1  3 

SCF  energy  of  C(  S)  and  C(  D)  respectively.  Finally,  the  total  energy 
+  3  • 

of  the  CH  Z  state  was  found  to  be  insensitive  to  the  4f  exponent.  No 

changes  were  made  in  the  H  basis  functions  in  the  enlarged  basis  set. 

All  of  the  reference  configurations  listed  in  Table  2  have  the  2a 

+  2  2  1 

orbital  doubly  occupied  although  one  of  the  asymptotes  C  (Is  2s2p  D)  + 


H  under  consideration  has  the  carbon  2s  orbital  singly  occupied.  In  the 
extended  calculations,  the  2a  orbital  was  included  in  the  valence  space 
and  the  configuration  list  was  not  limited  by  reference  to  a  set  of 
specific  configurations.  For  the  calculations  on  the  ^Z+  symmetry,  the 
valence  set  consisted  of  the  2a,  3a,  4 a  and  lrr  orbitals,  the  Rydberg  set 
included  the  5a,  6a,  7 a  and  2tt  orbitals  and  the  remainder  of  the  orbitals 
comprised  the  external  set.  The  configuration  list  included  all  config¬ 
urations  generated  by  distributing  (i)  4  electrons  in  the  valence  set  and 
(ii)  3  electrons  in  the  valence  set  and  one  in  the  Rydberg  set,  plus  all 
single  and  double  excitations  with  respect  to  the  configurations  in  (i) 
and  (ii).  For  calculations  on  the  ^2  symmetry,  the  Rydberg  set  consisted 
of  only  the  2tt  orbital.  The  la  core  orbital  was  kept  fully  occupied  in 
all  configurations.  This  procedure  resulted  in  15454  configurations  for 
the  V  symmetry  and  13063  configurations  for  the  ^2  symmetry. 

The  molecular  orbitals  used  in  the  extended  calculations  were 


determined  by  the  same  four  step  procedure  previously  described 


Calculations  of  the  five  lowest  £  states  and  four  lowest  il  states  were 
carried  out  at  a  limited  number  of  intemuclear  separations.  Results  of 
these  test  calculations  are  presented  and  compared  with  the  results  of  the 
calculation  described  previously. 

The  MCSCF  calculations  were  performed  using  a  program  developed  by 
J.  Hinze.  All  other  calculations  were  performed  using  the  ALCHEMY 
system  of  programs  developed  by  P.  S.  Bagus,  B.  Liu,  A.  0.  McLean  and 
M.  Yoshimine. 

Ill •  Results  and  Discussion 

A.  Potential  curves  and  asymptotes 

1  + 

In  Figure  2,  potential  curves  for  the  five  lowest  £  and  four  lowest 
states  of  CH+  are  presented.  The  energies  of  these  states  as  a  function 
of  intemuclear  distance  are  listed  in  Table  3.  Of  these  nine  states, 
only  the  X^E*",  aS  and  2^£+  states  are  bound  with  respect  to  their  asymp¬ 
totic  limits.  Dissociation  energies,  equilibrium  separations  and  spectro¬ 
scopic  constants  for  the  X  and  A  states  are  compared  in  Table  4  with  the  pre 
vious  Cl  calculations  of  Green  et  al.7  and  with  experiment.  Satisfactory 
agreement  with  both  sets  of  values  is  observed.  Discrepancies  for  the 
A^£  state  are  somewhat  larger  than  those  for  the  X*£+  state. 

Excitation  energies  of  the  four  asymptotes  considered  in  this  work, 

+  2  2 

relative  to  the  lowest  asymptote  C  (  P)  +  H(  S)  at  the  separated  atom 

limit  of  R  *  20a  are  given  in  Table  5.  The  excitation  energies  of  the 
o 


C(^D)  +  H  and  the  C(^S)  +  H  asymptotes  are  each  larger  by  0.4  eV  chan 
Che  experimental  values.  However,  the  calculated  splitting  between  the 
two  states  of  neutral  carbon  C(^D)  and  C(^S)  is  exactly  equal  to  the 
spectroscopic  value.  In  effect,  the  error  in  the  ionization  potential 
of  the  carbon  atom  is  -'■0.4  eV  in  this  calculation.  Excitation  energies 
are  given  critical  attention  in  che  comparison  between  these  values  and 
the  results  of  the  extended  Cl  calculations  presented  in  the  next  section. 

B.  Extended  Cl  Results 

The  extended  Cl  calculations  were  carried  out  at  the  R  ■  20a 

o 

asymptote  and  at  three  points  in  the  interaction  region.  Excitation 
energies  for  these  calculations  are  also  listed  in  Table  5.  Discrepancies 
in  excitation  energies  for  the  first  two  excited  limits  are  lowered  from 
0.4  eV  to  0.2  -  0.25  eV .  There  is,  however,  no  change  in  the  excitation 
energy  of  che  carbon  Rydoerg  limit,  C(^P)  +  H+.  Asymptotically,  the 
total  energy  from  the  extended  calculations  is  lower  by  3  millihartrees 
for  the  X  Z  state  and  by  2  mh  for  the  A  H  state. 

Interaction  potential  curves  were  compared  at  internuclear  distances 

of  1.9,  2.2,  and  3.4  a  .  For  all  nine  states,  the  maximum  energy  dis- 

o 

crepancy  observed  was  0.3  eV ;  all  others  are  £  0.2  eV.  In  the  extended 
calculations,  near  the  potential  minimum,  che  A^H  state  was  lowered  by 
-~0.2  eV  while  the  X^Z+  state  was  virtually  unchanged.  That  the  A  *11  is 
less  well  determined  in  the  main  calculations  reported  here  than  the 
X*Z+  state  is  further  confirmed  by  the  greater  discrepancy  with  respect 
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Co  experimental  values  of  the  spectroscopic  constants  in  Table  4  for  the 
Xl2  state  than  for  the  X^£  state.  However,  in  general,  energy  shifts  of 
this  magnitude  have  little  effect  on  the  qualitative  features  of  the 
potential  curves.  Thus  it  was  concluded  that  the  main  calculations 
reported  here  realistically  characterize  the  excited  states  of  CH+. 


C.  Discussion  of  Potential  Curves 

In  using  potential  curves  shown  in  Figure  2  to  interpret  physical 
processes,  the  most  reliable  information  is  obtained  by  shifting  the 
curves  so  that  their  asymptotes  have  the  spectroscopic  separation,  given 

in  Table  5. 


1  ■+*  1  •+• 

The  avoided  crossing  of  the  2  £  and  3  £  states  gives  rise  to  a 

humo  in  the  21£+  state  with  a  height  of  <0.5  eV  between  3  and  4a.  A 

o 


cursory  search  for  quasi-bound  levels  on  an  energy  grid  of  5  x  10 
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har trees,  led  to  three  resonances  at  0.296,  0.466  and  0.588  eV  above  the 

,  1  +  14 

separated  atom  limit.  The  2  £  state  appears  to  experience  a  small  /R 

attraction  at  large  R  with  an  energy  minimum  of  about  0.2  eV  at  R  »■  6.2  a 

This,  however,  is  well  outside  the  Franck-Condon  region  of  the  ground 

state . 

The  two  highest  ^£+  and  states  calculated  also  appear  to  interact 
The  4^£+  and  5^£+  states  (as  well  as  the  3^1  and  4^11  states)  are  very 


close  in  energy  at  R  «  3.5  a  ,  creating  ripples  in  the  potential  curves. 

o 

The  5^£+  also  appears  to  have  a  quasi-bound  vibrational  level  about 


2.8  eV  above  the  separated  atom  limit. 


It  is  clear  chat,  with  the  exception  of  the  quasi-bound  vibrational 

levels  mentioned  above,  no  new  band  spectra  of  CH  are  to  be  found.  The 

excited  states  will  be  important,  however,  for  photodissociation.  With 

the  asymptotes  shifted  to  agree  with  experiment,  the  vertical  excitation 

energies  at  R  =  2.1  a  from  the  v=0  level  of  the  ground  state  to  the 

o 

1  4*  1  +  i 

2  Z  ,  3  I  and  2  “  states  are  8.02,  12.92  and  13.58  eV,  respectively. 

These  s cates  may  be  significant  in  photodissociating  CH+  in  the  inter- 

+ 

stellar  medium.  In  comets  where  a  large  amount  of  CH  can  be  observed 
and  where  there  is  no  Lyman  cut-off  of  the  solar  photons,  even  higher 
states  may  be  important  in  the  photodestruction  of  CH+. 


D.  Transition  Moments 

The  electric  dipole  moments  for  transitions  between  the  ground  state 
and  the  2^2  ,  3^2+,  A  ^3  and  2^2  states  of  CH+  are  given  as  a  function  of 
R  in  Table  6  and  Figure  3.  Properties,  such  as  the  transition  moment, 
tend  to  be  more  sensitive  than  the  energy  to  the  quality  of  the  wavefunction 
used  to  calculate  them.  When  the  transition  moments  from  the  main  calcu¬ 
lations  were  compared  with  the  results  of  the  extended  Cl  calculations 
there  was  excellent  agreement  for  transitions  to  the  2^2+,  3^2+,  A^2  and 
2^1  states  but  fairly  large  discrepancies  for  transitions  to  the  higher 
states,  41'2+,  5^2+>  3^3  and  4^3.  Since  one  of  the  main  improvements  in 

the  extended  Cl  calculations  is  the  inclusion,  in  effect,  of  2s  excited 

+  2 

reference  configurations  and  since  the  C  (  D)  asymptote  to  which  the 
41I+  and  3^3  dissociate  is  a  2s  excited  configuration,  this  observation 


is  reasonable.  We,  therefore,  present  transition  moments  only  for  the 
lower  states,  which  include  all  states  assessible  in  the  interstellar 
medium. 

Our  values  for  the  X-A  transition  moment  are  generally  within  207.  of 

5 

those  of  Yoshimine  et  al."  which  was  used  to  predict  an  oscillator  strength 


for  transitions  between  the  lowest  vibrational  levels  of  each  state  in 

6 

excellent  agreement  with  a  subsequent  experimental  measurement.  At 

R  *  20  a  the  transition  moments  between  the  ground  state  and  the  first 
o 

two  excited  states  of  each  symmetry  are  zero,  consistent  with  the  fact  that 

+  2  2 

there  is  no  dipole  transition  possible  between  C  (  ?)  +  Hi,  S)  and  either 

C(^D)  +  H  or  C ( ''"S )  +  R+.  At  R  =*  2.1  a  ,  transition  moments  to  the  3^E+ 

o 

and  2^2  are  relatively  large,  in  agreement  with  the  finding  of  Watson 
9 

et  al.  The  sum  of  the  electronic  oscillator  strengths  to  all  the  excited 
states  presented  here  is  1.03.  There  are  many  peaks  and  dips  in  the  trans¬ 
ition  moments  for  the  higher  states  which  are  indicative  of  the  changes 
of  character  in  the  wavefunction  occurring  due  to  valence -Rydberg  mixing. 

Photodissociation  cross  sections  to  the  2^£+,  3^£+  and  2^2  and 
excited  states  will  be  reported  separate ly. ^  The  photodissociation  cross 


section  and  rates  for  the  A-X  transition  have  been  calculated  by  Uzer  and 

12  5 

Dalgarno,  using  the  transition  moments  of  Yoshimine  et  al. 
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Slater  Basis  Set  for  CH+ 
Extended  Calculations 
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Reference  Configurations  for  CT  Calculations 


Table  3 

POTENTIAL  CURVES  FOR  ‘l*  AND  'll  STATES  OF  CH+ 

GIVEN  IN  HARTREES  RELATIVE  TO  THE  ASYMPTOTIC  ENERGY  OF  EACH  STATE* 


x's* 

4’S+ 

5’r*- 

1.5 

-0.03335 

0.16712 

0.32063 

0.26491 

0.29191 

1.6 

-0.07786 

0.11777 

0.26916 

0.21644 

0 . 23624 

1.7 

-0.10867 

0.08192 

0.23058 

0.18132 

0 . 19408 

1.8 

-0.12924 

0.05614 

0.2012S 

0.15576 

0.16234 

1.9 

-0.14210 

0.03790 

0.17845 

0.13661 

0.13909 

2.0 

-0.14920 

0 .02530 

0.16021 

0.12126 

0.  12323 

2.1 

-0.15197 

0.01690 

0.14512 

0. 10849 

0.  11326 

2.2 

-0.15154 

0.01163 

0. 13221 

0.09807 

0.10726 

2.4 

-0. 14424 

0.00737 

0.11065 

0.08353 

0.10193 

2.6 

-0.13200 

0.00798 

0.09286 

0.07490 

0.10080 

2.3 

-0.11762 

0.01079 

0.07771 

0.06945 

0 . 10040 

3.0 

-0.10275 

0.01418 

0 . 06447 

0.06523 

0.09274 

3.2 

•0.08836 

0.01722 

0.05257 

0 . 06080 

0.07332 

3.4 

-0.07500 

0.01935 

0.04169 

0.05478 

0.0574O 

3.6 

-0.06293 

0.02027 

0.03183 

0.04600 

0.04781 

3.8 

-0.05225 

0.01979 

0.02331 

0.03623 

0.04273 

4.0 

-0.04297 

0.01788 

0.01644 

0.02753 

0.03996 

4.3 

-0.03151 

0.01290 

0.00934 

0.01738 

0 . 03820 

4.6 

-0.02270 

0.00692 

0.00S2O 

0.01039 

0.03748 

S.O 

-0.01433 

-0.00003 

0.00217 

0 . 00463 

0.03519 

5 .  5 

-0.00793 

-0.00541 

0.00020 

0.00117 

0.02858 

6.0 

-0.00441 

-0.00739 

-0.00083 

-0 . 00024 

0.01954 

7.0 

-0.00154 

-0.00711 

-0.00148 

-0.00075 

0.00168 

8.0 

-0.00068 

-O.OOSOl 

-0.00129 

-0.00063 

-0.01080 

20.0 

-0.00000 

0.14672 

0.19880 

0.35102 

0.36497 

R(*tj) 

ATn 

zfn 

3'n 

4’n 

1.5 

0 . 10645 

0.44383 

0.26283 

0.29568 

1.6 

0 . 06008 

0.38467 

0.20883 

0.24415 

1.7 

0.02710 

0.33672 

0.17097 

0.20622 

1.8 

0.00402 

0.29822 

0.14410 

0.17786 

1.9 

-0.01175 

0.26714 

0.12496 

0.15373 

2.0 

-0.02215 

0.24156 

0.11132 

0.12767 

2.1 

-0.02865 

0.21984 

0.10152 

0.10336 

2.2 

-0.03234 

0.20060 

0.09394 

0.08342 

2.4 

-0.03438 

0.16603 

0.07421 

0.06593 

2.6 

-0.03266 

0.13472 

0.05958 

0.06166 

2.8 

-0.02965 

0. 10713 

0.05349 

0.05991 

3.0 

-0.02655 

0.08373 

0.05126 

0.05841 

3.2 

-0.02267 

0.06644 

0.05273 

0.05635 

3.4 

-0.01992 

0.05146 

0.04998 

0.05425 

3.6 

-0.01746 

0.03954 

0.04393 

0.052S6 

3.8 

-0.01525 

0.03002 

0.03592 

0.05071 

4.0 

-0.01326 

0.02237 

0.02787 

0.04824 

4.3 

-0.01068 

0.01355 

0.01780 

0.04347 

4.6 

-0.00854 

0.00719 

0.01065 

0.03809 

5.0 

-0.00626 

0.00157 

0.00477 

0.03103 

5.5 

-0.00420 

-0.00223 

0.00115 

0.02335 

6.0 

-0.00281 

-0.00379 

-0.00031 

0.01710 

7.0 

-0.00132 

-0.00384 

-0.00083 

0.00708 

8.0 

-0 . 00068 

-0.00282 

-0.00063 

-0.00149 

20.0 

-0.00069 

0.14817 

0.35221 

0.37578 

*The  aaymptotic  energy  for  each  rtata  it  taken  to  be  the  calculated  energy  at  ft  ■  20  a0.  The 
tnergiei  for  ft  *  20  ea  ire  given  relative  to  the  calculated  x'S*  atymotote  of  -37.3703889? 
hanreet.  1  aQ  •  0.52918  A;  1  hartree  ■  27.2116  ev. 
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Table  4 


Spectroscopic 

+  1  + 

Constants  for  CH  X  Z 

and  A1!!  States 

This  Work 

Green  et  al. 

Experiment 

(Ref.  7) 

(Ref.  13) 

x1! 

R  (a  )  2.133 

e  o 

2.136 

2.137 

D  (eV)  4.140 

e 

4.11 

4.27a  =0.02 

AS. (cm‘1)2743 
h 

2740 

AlZ 

R  (a  )  2.381 

e  o 

2.332 

2.333 

D  (eV)  0.936 

e 

1.069 

1.29° 

(cm  )  1526 
e 

1865d 

T  25690 

e 

24970 

24033® 

a)  D  °  +  uj  /2  using  values  in  Ref.  13 
o  e 


b)  Using  4  vibrational  levels 

c)  Using  experimental  T  ,  D  0  for  X  state  and  uu  /2  for  A  state 

oo  e 

d)  Applies  to  v  s  2 

e)  Using  experimental  T  and  a  /2  for  X  and  A  states. 

o  e 
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H*  2  2  & 

Excitation  Energies  Relative  to  C  (  P)  +  H(  S)  in  eV 


V 

Extended 

Calculations 

'a 

Extended 

Calculations 

Spectroscopic 

1  + 

C(  P)  +  H 

9.93 

9.93 

10.22 

10.22 

10.023 

+  2  2 

C  (  D)  +  H(  S) 

9.55 

9.44 

9.58 

9.44 

9.290 

1  + 

C(  S)  +  H 

5.41 

5.28 

5.022 

1  + 

C(  Q)  +  H 

3.99 

3.82 

4.03 

3.83 

3.602 

a)  Asymptote 

(R»20a  ) 

for  X^£+  state 

• 
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Table  6 

CH+  Transition  Dipole  Moments 


R 
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i 

[  ; 
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0.6731 
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-0.5588 

-0.3611 

0.6873 

i  i 

1.8 

0.2165 

-0.6563 

-0.3373 

0.6873 

i 

1.9 

0.1965 

-0.7592 

-0.3123 

0.6781 

" 

1 

2.0 

0.1747 

-0.8653 

-0.2881 

0.6592 

;  1 

2.1 

0.1507 

-0.9712 

-0.2634 

0 . 6287 

2.2 

0.1244 

-1.0738 

-0.2389 

0.5862 

1 

2.4 

0.0633 

-1.2591 

-0.1919 

0.4762 

2.6 

-0.0110 

-1.4105 

-0.1492 

0. 3653 

2.3 

-0.1012 

-1.5277 

-0.1128 

0.2782 

. 

3.0 

-0.2102 

-1.6118 

-0.0833 

0.2116 

3.2 

-0.3410 

-1.6617 

-0.0563 

0.1561 

3.4 

-0.4958 

-1.6736 

-0.0380 

0.1137 

i 

i 

3.6 

-0.6744 

-1.6415 

-0.0246 

0.0790 

3.8 

-0.8694 

-1.5598 

-0.0151 

0.0507 

4.0 

-1.0600 

-1.4308 

-0.0088 

0.0276 

4.3 

-1.2777 

-1.1888 

-0.0032 

0.0013 

\  4-s  ' 

5.0 

-1.3743 

-0.9608 

-0.0008 

0.0165 

-1.3510 

-0.7363 

-0.0003 

0.0301 

| 

j  »  5.5 

-1.1880 

-0.5576 

-0.0008 

0.0359 

1 

s 

6 . 0 

-0.9721 

-0.4409 

-0.0015 

0.0346 

i 

7.0 

-0.5783 

-0.2369 

-0.0017 

0.024S 

i 

3.0 

-0.3173 

-0.1820 

-0.0012 

0.0151 

i 

20.0 

0.0002 

0.0003 

0.0000 

0.O0C1 
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Figure  Captions 


1  +  1 

Figure  1:  Separated -a too  limits  giving  rise  to  £  or  ^II  states 
of  CH+.  Experimental  energy  separations  are  shown  as 
taken  from  tables  compiled  by  C.  E.  Moore,  MBS 
circular  467,  1949. 

Figure  2:  Potential  curves  of  CH+ :  five  of  ^E+  symmetry  and 
four  of  Ln  symmetry. 

Figure  3:  Transition  moments  from  the  ground  state  to  the 

2^E+,  3^E+,  AlH  and  2^H  states  of  CH+  as  a  function 
of  intemuclear  separation.  The  relative  signs  of 
the  transition  moments  are  arbitrary. 
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ABSTRACT 


A  simple  model  that  has  been  used  extensively  by  Los  and  coworkers 
to  treat  ion  pair  formation  in  collisions  of  alkali  atoms  with  diatomic 
molecules  is  extended  to  include  continuum  coupling  via  a  competing 
Penning  ionization  channel.  This  extended  model  is  then  used  to  calcu¬ 
late  the  differential  cross  sections  for  ion  pair  formation  for  the 
* 

system  Ar  +  I^  over  the  energy  range  28-154  eV  and  to  compare  with  a 

previous  treatment  of  K  +  I2 •  In  the  absence  of  significant  competition 

from  continuum  processes,Ar  is  expected  to  behave  in  a  manner  similar 

to  K,  since  the  active  electron  is  an  unpaired  4s  electron  in  both  cases. 

* 

We  perform  model  calculations  for  Ar  +  I^  to  investigate  the  effects  of 

varying  the  potential  curves  and  charge  exchange  matrix  elements  and  of 

including  a  continuum  coupling  function  F  (R) .  Comparison  with  previous 

* 

calculations  for  K  +  I^  suggests  increased  repulsion  on  the  Ar  -  I^ 
surfaces  relative  to  those  of  K  -  I  .  The  competing  mechanisms  of 
excitation  transfer  and  Penning  ionization  may  have  a  small  effect  upon 
the  ion  pair  angular  distributions. 
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I  INTRODUCTION 

The  study  of  a  pair  of  generally  similar  reactions  enables  one  to 
isolate  subtle  effects  that  arise  specifically  from  the  differences 
between  the  two  systems.  Recent  measurements  of  the  differential  cross 
section  for  ion  pair  formation  in  collisions  between  halogen  molecules  and 
alkali  atoms  or  metastable  rare  gas  atoms  provide  such  an  opportunity. 

Over  the  energy  range  30-150  eV,  all  of  the  major  features  observed^- 
in  the  angular  distribution  for  the  reaction 

Ar*  (3p2,0}  +  l2  -  Ar+  +  h  (1) 

2-4 

can  be  correlated  with  those  observed  in  the  analogous  reaction 

K  +  I2  -  K+  +  I2*  .  (2) 

However,  there  are  differences  in  the  ratios  of  the  intensities  of  various 

features  and  in  the  position  of  the  rainbow  angle.  The  general  similarity 

*  3 

is  consistent  with  the  fact  that  both  K  and  Ar  (  P)  have  an  active, 

unpaired  4s  electron.  The  observed  differences  between  the  two  reactions 

* 

should  be  related  to  the  inner  shell  vacancy  in  Ar  and  to  the  additional 

competing  channels  energetically  accessible  to  the  system,  including 

Penning  ionization: 

*  + 

Ar  +  I2  -  Ar  +  I„  +  e  .  (3) 

In  this  paper  we  present  model  calculations  designed  to  elucidate  the 
causes  of  the  observed  differences  in  the  scattering  data. 

C-2 
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Ion-pair  formation  in  the  alkali-halogen  system  has  been  extensively 

investigated"5  ^  both  experimentally  and  theoretically.  These  reactions 

5,6 

have  been  modelled  quite  successfully  by  Los  and  coworkers  using 

8  9 

classical  trajectory  surface-hopping  techniques  ’  that  demonstrate  the 

essential  features  of  the  reaction  mechanism  leading  to  ion-pair  production. 

Here  we  have  modified  and  extended  those  model  calculations  to  include 

coupling  to  the  continuum  (reaction  (3)).  We  estimate  the  autoionization 

rate  T  (R)  using  a  simplified  form  of  a  technique  proposed  recently  by 
10 

Miller  and  Morgner.  We  are  therefore  able  to  perform  quantitative 


calculations  that  examine  the  effect  of  competing  channels  upon  the  ion 
pair  distributions  in  reaction  (1) . 

In  Section  II  we  present  the  details  of  the  model  calculations  and 
discuss  the  inclusion  of  continuum  coupling.  Our  results  are  presented 
and  discussed  in  Section  III,  and  Section  IV  contains  a  summary. 


i 
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I  (  , 


II  THEORY 


A.  Summary  of  Che  Model 

Me  have  adopted  Che  theoretical  approach  developed  by  Los  and  co- 

2  3  6  3  9 

workers,  ’  ’  which  is  a  classical-trajectory ,  surface-hopping  model  ’  ’ 

j 

with  the  following  additional  approximations:  ' 

1.  Simple  analytic  potentials  are  constructed  from  pairwise 
interactions  for  the  covalent  electronic  state  K  +  I^  and 
for  the  ionic  state  K.  +  I^  • 

2.  The  cross  sections  are  obtained  as  the  average  of  cross 
sections  for  various  orientations  of  the  target  molecule. 

The  axis  of  this  target  is  assumed  not  to  rotate  during 
the  collision.  Molecular  vibration,  however,  is  included, 
as  explained  below. 

3.  For  each  trajectory,  the  classical  vibration  of  the  molecule 
is  calculated  numerically  assuming  Morse  potentials,  and 

the  deflection  of  the  projectile  is  obtained  analytically 
using  classical  perturbation  theory. 

An  essential  feature  of  the  model  is  that  the  ionic  and  covalent 
* 

surfaces  cross  at  an  Ar  -I^  distance  R  that  depends  on  the  vibrational 
coordinate  of  Che  molecule.  This  fact  influences  the  collision  in  the 
following  way.  The  system  begins  asymptotically  on  che  covalent  surface. 

When  the  particles  reach  Che  location  of  the  crossing  with  che  ionic 
surface,  the  Landau-Zener  formula^  is  used  to  compute  the  probability  of 
switching  to  the  other  surface.  As  shown  in  Figure  1,  there  are  two 
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possible  paths  to  ion-pair  formation.  The  surface  crossing  may  occur  on 

the  inward  part  of  the  trajectory  (the  dashed  line)  or  on  the  outward  part 

(shown  by  the  dotted  line) .  A  transfer  to  the  ionic  surface  will  initiate 

an  expansion  of  the  bond  since  has  a  larger  equilibrium  intemuclear 

distance  than  I?.  If  this  transfer  occurs  on  the  incoming  trajectory,  the 

resulting  I?  expansion  will  yield  an  effective  vertical  electron  affinity 

that  increases  with  time.  This  will  increase  the  radius  R  at  which  the 

second  surface  crossing  is  encountered  on  the  outward  part  of  the  trajectory. 

Since  the  matrix  element  connecting  the  two  surfaces  decreases  strongly 

with  increasing  R,  the  probability  of  a  diabatic  path  through  the  second 

crossing  may  be  significantly  larger  than  at  the  first  crossing.  These 

1-6 

effects  have  already  been  thoroughly  discussed  in  the  literature. 

The  analytic  form  of  the  potential  surface  we  have  used  differs 
somewhat  from  the  form  used  in  reference  2.  Ours  is  defined  in  terms  of 
Figure  2  to  be 

voo  ■  V(V  +  V(V  +  Vr)  (4) 


V<V  -  ^  *  V(R2)  -  2**  ♦  VD 


Vg(r)  and  v^(r)  are  Morse  potentials  for  I  and  respectively.  The 
same  Morse  parameters  were  used  as  in  reference  2,  namely,  0  *  1.54  e^, 

3  -  1.87  A-1,  and  r^  -  2.67  A  for  I  ,  and  D'  -  1.02  eV,  3'  «  1.23  JTl, 
and  zj  ■  3.20  A  for  .  V(R  )  is  a  repulsive  term  of  the  Born-Mayer  form, 

-BR 

V(*t)  -  A  e  i  (6) 
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where,  foLlowing  reference  2,  we  initially  use  A  ■  6.4  x  10  eV  and 
•  -1 

3  ■  4.762  A  .  Note  that  for  simplicity  the  same  constants  A  and  B 

are  used  here  for  the  ionic  and  covalent  terms.  The  potentials  we 

use  are  much  simpler  than  those  of  reference  2.  Since  ab  initio 

calculations  are  not  available,  we  choose  to  use  potentials  with  as 

few  adjustable  parameters  as  possible. 

The  potential  surfaces  defined  by  Eqs.  (4)  and  (5)  do  not  reflect 

★ 

the  multiplicity  of  states  arising  from  the  inner  shell  vacancy  in  Ar  . 
Our  simplified  model  is  equivalent  to  assuming  chat  all  of  the  covalent 
surfaces  are  degenerate,  and  that  all  of  the  ionic  surfaces  are 
degenerate.  This  assumption  appears  reasonable  at  large  values  of  R, 
where  the  surfaces  cross.  In  this  region  we  expect  no  inner  shell  re¬ 
arrangement  processes.  The  core  electrons  should  have  only  a  minor 
influence  on  the  potential  shapes  and  coupling  matrix  elements.  This 
fact  is  indeed  the  major  reason  for  the  strong  analogy  to  alkali 

reactions  and  the  justification  for  our  use  of  the  single  surface 

13 

crossing  and  Che  coupling  matrix  element  angular  dependence  (cos# 

in  Eq.  (7)  below)  that  are  appropriate  to  alkali  systems.  The  inner 

■k 

shell  vacancy  in  Ar  would,  however,  be  expected  to  influence  surface 
parameters  at  much  smaller  R  values  and  thereby  could  affect  the  shape 
of  the  differential  cross  sections.  We  will  return  to  this  point  later. 


We  use  the  same  form  for  the  coupling  matrix  element  as  in  reference  2 

* 

.  ^  -c^R 

v12^R,r'3'1  *  c]_^1  *  EA(r)]^  Re"  cos  a  (7) 

where  R*  (1^  +  EA(r)^).  (8) 

y  z 

I  is  the  ionization  potential  of  the  projectile  and  EA(r)  is  the  vertical 
electron  affinity  of  the  molecule  as  a  function  of  the  separation  r. 

o 

in  eV;  I  and  EA  in  a.u.;  R  in  A;  c^  in  A  ).  Note  that  the  diagonal 

terms  and  V  are  orientation-dependent  because  they  are  defined  in 

terms  of  R  and  R  .  Note  also  that  the  angular  dependence  of  V  forces 
12  12 

13 

it  to  be  zero  at  the  angles  required  by  symmetry.  Initially  the 

i-l  2 

constants  c^  *  0.28  eV  and  c^  ■  0.65  A  were  used  for  both  reaction  (1) 
and  reaction  (2). 

For  each  orientation  of  the  target  (in  a  space-fixed  frame),  we  wish 
to  calculate  the  two  branches  of  the  classical  deflection  function. 

These  two  branches  give  the  scattering  angle  as  a  function  of  impact 
parameter  for  trajectories  in  which  the  surface  crossing  occurs  on  the 
inward  or  outward  part  of  the  trajectory.  (If  the  crossing  occurs 
twice,  or  not  at  all,  ion-pair  formation  does  not  occur  and  the  trajec¬ 
tory  is  not  counted.)  For  each  impact  parameter,  the  scattering  angles 
are  therefore  calculated  for  two  sequences  of  events.  In  both  cases, 
the  1^  molecule  begins  at  rest  on  the  neutral  potential  curve  Vg(r)  at 
the  equilibrium  distance  r@,  and  the  projectile  follows  a  straight-line 
path  with  constant  velocity  throughout  the  collision.  For  the  ionic 


branch  of  Che  deflection  function,  Che  system  is  assumed  to  switch  to 

the  ionic  potential  at  the  first  crossing.  Then  the  1^  molecules 

begin  to  move  classically  on  the  potential  v^(r),  beginning  at  rest  at 

r  *  r  (I  ) .  The  location  of  the  second  crossing  is  determined 
e  2 

numerically  by  monitoring  the  oscillatory  motion  of  the  target  and  the 
rectilinear  motion  of  the  projectile.  For  the  covalent  branch  of  the 
deflection  function,  the  molecules  do  not  begin  to  vibrate  until  the 
crossing  is  encountered  a  second  time,  so  that  the  first  and  second 
crossings  occur  at  the  same  intermo lecu lar  separation  R.  For  both 
branches  of  the  deflection  function,  the  orientation  angle  y 
is  computed  at  both  crossings.  The  orientation  influences  the 
transition  probability  through  the  factor  cos  y  in  Eq.  (7). 

The  preceding  discussion  shows  how  the  appropriate  potential  is 
determined  at  each  point  along  a  given  trajectory.  The  scattering 
angle  is  then  determined  by  integrating  the  cumulative  momentum  trans¬ 
fer  perpendicular  to  the  rectilinear  motion,  according  to  standard 
14 

formulas  of  classical  perturbation  theory.  The  deflection  is  easily 
calculated  analytically  because  the  functions  in  equations  (4)  and  (5) 
are  written  as  sums  of  two-body  potentials. 

In  order  to  obtain  the  cross  section,  it  is  also  necessary  to 
calculate  the  probability  that  the  system  will  follow  a  given  trajectory 
This  is  dona  using  the  Landau-Zener  formula. ^  The  total  probability 


for  a  given  sequence  of  surface  crossings  or  avoided  crossings  is  the 


product  of  the  probabilities  for  each  event  separately.  This  is  easily 
calculated  once  the  locations  of  the  two  crossing  points  are  known.  The 
result  for  the  cross  section  (for  a  particular  orientation  of  the  molecule) 
can  be  summarized  in  the  following  formula: 


dc  . 

“  sin  9 
dO 


r- 

j 


ion  . 


|d0lon/db 


j' 


/PC°V(bj  Q 
ldSC0V/db1 


(9) 


ion  cov 

P  (P  )  is  the  overall  probability  of  following  an  ionic  (covalent) 
trajectory  for  a  particular  impact  parameter  b.  The  sums  are  understood 
to  be  over  all  impact  parameters  that  lead  to  the  scattering  angles  9  or 
-9.  Although  eq.  (9)  presents  the  formal  definition  of  the  cross  section, 
we  found  it  useful  in  practice  to  use  a  conventional  histogram  procedure. 


B.  The  Penning  Ionization  Channel. 

1.  Estimating  P(R)  for  molecular  systems. 

F(R) A  is  the  rate  at  which  particles  are  lost  from  an  initial 
channel  into  the  ionization  channel  via  the  reaction  (3).  P  is  given  by^’^ 


r  -  2rrp | H 


(10) 

(11) 


where  H  is  the  Hamiltonian,  }  ,  is  the  initial  bound  state  wave  function 

initial 

*  + 

(Ar  +  I«) ,  i..  ,  is  the  continuum  configuration  (Ar  +  I„  +  e  )  and  p 

l  tinai  z 


is  Che  density  of  continuum  states.  If  the  continuum  wave  function  is 
normalized  to  unit  amplitude  asymptotically,  then 

2np  -  4/k  (12) 


where  the  kinetic  energy  of  the  ejected  electron  is,  in  atomic  units, 


E  -  1/2  k  . 


10 


Miller  and  Morgner  have  suggested  that  the  matrix  element  can 

be  written  as  the  product  of  an  overlap  factor  and  a  charge  exchange 

* 

factor.  Specializing  to  our  case,  Ar  +  I^,  this  would  be  the  product 


H  as  S  x  H 
fi  fi  fi 


(o) 


(13) 


where  S  is  the  overlap  between  the  valence  electron  orbital  on  Ar  and 

(o) 

the  final  state  continuum  orbital,  and  R  is  the  matrix  element  for 


the  one-electron  process 

Ar+  +  I  —  Ar  +  I.+ 
2  2 


(14) 


Miller  and  Morgner10  pointed  out  that  can  estimated  using  the 


semi-empirical  correlation  formula  of  Olson  et  al.: 
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[Hfi(o)l  as  OR  e"R/C 


(15) 


where  R  is  the  internuclear  separation  in  a.u.,  and 


1/0  ■  h  [1J2  ♦  i1 


(16) 


G  -  (1I2  V  /C 


(17) 


where  IT  and  I  are  the  ionization  potentials  of  I,  and  Ar  in  a.u.. 
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respectively.  |h  [  is  then  given  also  in  a. 


19 

Based  on  past  experience  calculating  continuum  overlap  factors 
^fi*  We  ^ave  c^e  Miller-Morgner  formula  by  setting  S  «  1. 

This  approximation  enables  us  to  evaluate  H  without  numerical  integra¬ 
tion.  Furthermore,  since  the  potential  surfaces  are  not  in  general  known, 
we  calculate  the  density  of  states  2rrp  from  the  energy  of  the  ejected 

electron  as  R  -  ®.  This  information  is  available  from  spectroscopic  data. 

10 

Miller  and  Morgner  estimate  that  their  formula  is  reliable  within 

a  factor  of  three  to  five.  We  expect  that  our  additional  approximations 

will  degrade  the  accuracy  somewhat,  but  the  result  should  certainly  still 

be  a  reasonable  order-of-magnitude  estimate. 

2.  Ionization  Probability  as  a  Function  of  Impact  Parameter. 

The  autoionization  rate  obtained  in  the  previous  section  has 
2  -Ml 

the  form  ?(R)  *  R  e  .  However,  we  have  found  that  over  the  important 
range  of  R  probed  in  the  collisions,  the  numerical  values  of  T  can  be 
well  fit  (—107.)  by  the  simpler  form 


r(R)  «  a1  e 


-a2R 


(18) 


This  is  the  form  of  T  that  has  been  assumed  in  many  previous  semi -empirical 
20»21 

studies. 

We  have  found  chat  for  a  T  of  the  form  of  Eq.  (18),  a  very  simple 
formula  can  be  obtained,  in  the  perturbation  limit,  for  the  ionization 
probability  as  a  function  of  impact  parameter.  Since 


C-ll 


che  translational  motion  is  already  treated  using  classical  perturbation 
theory,  a  perturbation  assumption  for  the  continuum  coupling  is  con¬ 
sistent  with  the  rest  of  the  model. 

Let  b  be  the  impact  parameter.  Then  the  survival  probability  P(b) 
that  ionization  does  not  occur  on  a  trajectory  of  impact  parameter  b  is 


P(b)  -  exp[-2F(b)] 


(19) 


where  in  the  perturbation  limit 


F(b) 


rOOdR 


a  MSlaS 

I  *v/l-(b2/RZ) 


(20) 


and  v  is  the  incident  velocity, 
o 


Substituting  Eq.  (18)  into  Eq.  (20),  integrating  by  parts,  and 


substituting  x  ■  R/b,  we  obtain 

F(b) 


-a0bx 


ftv 


dx 


(21) 


The  integral  can  be  represented  in  terms  of  the  modified  Bessel  function 


f  * 
'  , 


.  Then 

F<b)  '  S'  h<hb)  ■ 

o 

Further  simplification  is  possible  because  of  the  relation 

\  -0.932  x 

x  K^(x)  *  2.076  e 


(22) 


(23) 


which  is  valid  to  about  107.  for  2  £  x  £  12.  The  final  result  is 
2.076a. 


F(b) 


2  0 


1  # -0.932  a2b 


(26) 


2 


Formulas  (19)  and  (24)  may  now  be  incorporated  into  the  model  in 
the  following  way.  For  each  orientation,  a  deflection  function  is  cal¬ 
culated  that  gives  the  angular  scattering  as  a  function  of  impact  parameter. 
When  the  cross  section  is  calculated  from  a  particular  branch  of  the 
deflection  function,  a  factor  is  first  included  to  represent  the  probabil¬ 
ity  that  the  particles  follow  the  correct  sequence  of  curve  crossings. 

The  result  is  then  multiplied  by  P(b),  which  gives  the  probability  that 
the  trajectory  is  completed  without  loss  to  ionization. 
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Ill  RESULTS  AND  DISCUSSION 


A.  K  +  I, 

2-4 

This  reaction  has  already  been  treated  in  detail.  Here  we 

demonstrate  that  the  use  of  the  potentials  defined  in  Section  II-A 

leads  also  to  satisfactory  agreement  with  experiment.  In  Figure  3, 

we  compare  the  calculated  and  experimental  values  at  a  lab  energy  (E^) 

of  60  eV.  The  agreement  is  not  surprising  since  our  potentials  are 

nearly  the  same  as  chose  fitted  to  the  data  in  reference  2.  The  major 

difference  is  that  we  have  neglected  the  small  van  der  Waals  term. 

* 

B.  Ar  +  1^:  General  Approach 

1  * 

Gillen  ec  al.  pointed  out  that  their  results  for  Ar  +  I?  were 

n 

similar  to  those  of  Aten  et  al.“  for  K  +  I  .  The  important  differ- 

* 

ences  are  that  the  covalent  peak  is  relatively  somewhat  smaller  for  Ar 
+  1^  than  for  K  +  I0  and  that  the  rainbow  in  the  ionic  peak  occurs  at 
smaller  angles.  The  first  difference  could  be  caused  by  the  continuum 
coupling  in  the  metastable  system  or  by  differences  in  the  potentials  and 
charge  exchange  matrix  elements.  The  second  difference  should  be  asso¬ 
ciated  with  a  difference  in  the  shape  of  the  potential  surfaces. 

In  this  and  the  fo-llowing  subsections  we  describe  our  quantitative 

calculations  to  investigate  possible  causes  for  the  experimental  dif- 

★ 

ferences.  With  only  a  qualitative  knowledge  of  the  Ar  +  I  potential 
surfaces,  we  do  not  expect  to  determina  uniquely  the  coupling  parameters 


★ 

in  this  system.  Hence,  our  approach  to  Ar  +  I,  ion-pair  formation 
calculations  is  to  start  with  the  K  +  1^  parameters  and  determine  the 
changes  in  the  distributions  caused  by  individually  altering  the  continuum 
width  (7-0  initially),  the  surface  coupling  parameters,  and  the  inter¬ 
action  potentials.  The  results  should  yield  insight  into  the  possible 

★ 

causes  of  the  differences  observed  between  Ar  +1,  and  K  +  I_  in  the 
ion-pair  channel. 

* 

We  initially  performed  trial  calculations  for  Ar  +  I,  assuming 

7  -  0  ,  and  using  the  same  potential  parameters  as  for  K  +  I,.  The 

crossing  point  differs  slightly  for  this  system  because  the  ionization 
★ 

potential  of  Ar  is  slightly  smaller  than  that  for  K.  This  change 

increases  somewhat,  and  thereby  decreases  V  slightly,  but  the 

effect  on  the  calculations  was  quite  small.  The  calculated  peak  ratios 

and  the  position  of  the  rainbow  angle  were  essentially  the  same  as  those 

calculated  for  K  +  I^,  and  these  values  did  not  agree  with  the  data  for 
* 

Ar  +  I2. 

★ 

C.  Ar  +  I^:  The  Continuum  Channel. 

We  then  examined  the  effect  of  including  continuum  coupling.  It 

should  be  noted  that  7(R)  is  assumed  to  affect  only  the  covalent  potential. 

22 

This  is  consistent  with  the  conclusions  of  Hultsch  et  al.  on  similar 
systems.  One  argues  that  Penning  ionization  of  the  covalent  electronic 
state  is  initiated  by  a  one-electron  process,  namely  the  "exchange" 
process,  in  which  an  electron  from  the  ground  state  atom  jumps  into  the 


./  i 


hole  of  the  excited  atom.  There  is  then  a  large  overlap  between  the 
valence  electron  of  the  excited  atom  and  the  continuum  orbital.  In 
contrast  for  a  system  on  the  Ar+  +  surface,  continuum  coupling 
necessitates  a  simultaneous  two-electron  rearrangement;  such  a  process 
would  be  expected  to  have  a  much  lower  probability. 

For  the  ground  electronic  state  of  I^+,  which  is  9.3  eV  above  I  ,  we 
follow  the  prescription  of  Section  II-B  and  obtain  (for  T  and  R  in  atomic 
units) 


r(R) 


1.752  R2 


-1.902R 

e 


(25) 


For  R  in  the  range  4  -  12  a  ,  this  is  well  fit  bv  the  form 

o 


-  __  ,,  , ,  -1.614R 

A  (R)  as  11.16  e  a.u. 

o 


(26) 


We  performed  calculations  using  this  form  of  T(R)  and  keeping  the 
other  potential  parameters  the  same.  The  result  is  that  the  continuum 
coupling  is  not  large  enough  to  cause  a  significant  effect  on  the  ion-pair 
distributions . 

The  charge  transfer  matrix  element  estimated  from  Eq.  (15)  is  a 


one-electron  matrix  element.  In  fact  I"(R)  should  be  larger  than  T  (R) 

o 

2 

to  account  for  the  two  equivalent  electrons  in  the  tt,  orbital. 

T,g 

2  2 

Moreover,  the  I  molecule  has  a  valence  configuration  (a  5p)  (rr  5p) 

2  S  5  f u 

2  2  2 

(it,  5p)  (tt  5p)  (rr,  5p)  ,  and  these  ten  electrons  have  vertical 
I> J  ■? » 8  T>  8 
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ionization  potentials  between  9.3  and  13.0  eV.  In  the  spirit  of  our 
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one-electron  estimate  of  T  ,  we  should  also  calculate  approximate  continuum 

o 

coupling  widths  for  the  other  electrons. 

Although  the  ionization  potential  for  electrons  in  the  lowest  two 

* 

of  these  orbitals  is  larger  than  the  Ar  excitation  energy.  Miller  and 
10 

Morgner  have  shown  that  the  electronic  energy  transfer  to  a  pseudo¬ 
continuum  Rydberg  state, 

Ar*  +  I2  -  Ar  +  I2**  (27) 

(core-excited  in  this  case)  can  be  treated  in  an  equivalent  way  to  a  true 

continuum  interaction.  The  details  of  the  formula  are  slightly  different; 

however,  the  calculated  coupling  widths  are  found  to  be  quite  similar  to 

those  for  states  whose  ionization  potentials  are  lower  than  the  excitation 

energy  of  the  incoming  state.  Assuming  the  contributions  for  all  ten 

electrons  can  be  added  (a  reasonable  first  approximation),  we  obtain  a 

resultant  T(R)  that  is  10-15  times  larger  than  V  (R)  for  R  values  between 

o 

3  and  11  a.u. 

The  effects  of  various  coupling  functions  T(R)  on  the  calculated 

angular  distribution  for  reaction  (1)  at  28  eV  are  shown  in  Figure  4. 

Continuum  effects  are  generally  more  pronounced  at  lower  energies,  and 

E  -  28  eV  is  the  lowest  energy  studied.  Clearly  T  -  10  T  ,  which 
o  o 

approximates  the  summed  ten  electron  continuum  coupling  estimated  above, 
yields  a  noticeable  depletion  of  the  covalent  peak.  However,  P  >  100  P 

o 

seems  required  before  the  intensity  ratios  of  the  covalent  and  ionic 


peaks  are  in  reasonable  agreement  with  experiment.  Since  our  approximate 


estimates  of  F  are  relatively  crude,  it  is  possible  that  F  is  indeed  an 

order  of  magnitude  larger  than  estimated.  Hence,  the  continuum  coupling 

effect  could  make  some  contribution  to  the  attentuation  of  the  covalent 

peak  in  reaction  (1).  However,  the  larger  T  values  (100  F  ,  1000  I*  ) 

o  o 

give  covalent  peaks  that  are  strongly  attenuated  at  low  T  values;  the 
angular  positions  of  the  calculated  peaks  then  disagree  considerably  with 
the  data. 

* 

D.  Ar  +  Sensitivity  to  the  Coupling  Matrix  Element  V^. 

* 

Another  possible  cause  for  the  smaller  covalent  peak  in  Ar  +  I, 

could  be  a  larger  coupling  matrix  element  V  We  have  performed  calculations 

in  which  F  *  0  and  the  constant  c^  of  V  is  varied  [see  Eq.  (7)]  .  We  find 

that  a  significant  change  in  V  from  the  value  for  K  +  is  needed  to 

★ 

cause  the  desired  change  in  Ar  +  peak  ratios.  Figure  5  shows  the 

experimental  angular  distribution  at  E  »  51  eV,  compared  with  three 

o 

calculations  assuming  V  multiplied  by  factors  of  1.0,  2.5,  and  5.0 

relative  to  K  +  I^.  Even  though  the  peak  ratios  are  accurately  matched 

by  scaling  V  ^  by  approximately  2.5,  we  feel  that  this  is  an  unreasonable 

* 

change  considering  the  expected  similarities  between  K  and  Ar  .  Although 
the  potential  surfaces  for  the  two  reactions  have  quite  different 
symmetries,  these  differences  are  associated  with  the  synsnetry  of  the 
Ar  core.  At  the  very  large  crossing  radius  (~5.7  X)  between  the  ionic 


8 


and  covalent  surfaces,  the  properties  of  the  core  should  play  a  minor  role 

* 

in  the  electron  transfer  process  and  the  outer  Ax  electron  should  behave 
like  a  4s  electron  in  K. 

* 

E.  Ar  +  1 2:  Sensitivity  to  V^. 

It  is  obvious  from  Figure  S  that  the  calculated  rainbow  angle  is  not 
in  good  agreement  with  the  experimental  data.  None  of  the  modifications 
detailed  above  has  altered  the  calculated  rainbow  position.  Since  its 
position  is  most  sensitive  to  the  veil  depth  in  the  ionic  channel,  we 
investigated  modifications  to  the  potential  parameters  necessary  for 
achieving  a  match  to  the  experimental  rainbow  angle.  In  our  model,  the 
well  depth  is  controlled  by  the  choice  of  the  parameters  A  and  B  in 
Eq.  (6).  We  found  that  varying  these  parameters  to  decrease  the  well 
depth  caused  the  rainbow  to  shift  to  smaller  angles,  as  expected.  Another 
effect  of  this  change  was  a  reduction  in  the  size  of  the  covalent  peak 
relative  to  the  ionic  peak.  Figure  6  shows  the  effect  of  changes  in 
the  potential  parameters  upon  the  ion-pair  distributions.  The  modification 
that  gave  a  reasonable  fit  to  the  ionic  rainbow  angle  corresponds  to  an 
Ar+  -  Ij  "well-depth"  that  is  -0.5  eV  shallower  than  that  of  K+  +  I  " 

(at  the  same  I-I  intemuclear  separation) . 

For  comparison,  there  are  good  ab  Initio  calculations  for  the 

24 

diatomic  potentials  in  ArF,  KrF,  and  the  Xa  halides.  For  each  of 
these  systems  the  well  depths  of  the  various  ion-pair  states  are 


between  0.3  eV  and  0.7  eV  shallower  than  those  of  the  analogous  alkali 

halide  molecule.  With  spin-orbit  coupling  included,  the  three  ion-pair 

states  calculated  for  each  rare-gas  halide  pair  differ  by  less  than 

—0.2  eV  in  well  depth  and  0.1  i  in  R  from  each  other.  If  a  similar 

result  obtained  for  the  triatomic  systems,  our  method  of  treating  the 
* 

Ar  +  I  system  by  a  single  (average)  ion-pair  surface  intersecting  a 

single  covalent  surface  would  be  a  reasonable  approximation  to  the 

average  scattering  on  the  relevant  surfaces  and  our  resultant  0.5  eV 
+  •  *+*  — 

difference  in  Ar  -  I?  and  K  -I^  well  depths  would  appear  to  be  reason¬ 
able. 

* 

F .  Ar  +  I  :  Cone lusions . 

2 

Our  calculations  clearly  indicate  that  the  ton-pair  angular  distri¬ 
butions  for  reaction  (1)  depend  quite  sensitively  on  the  shapes  of  the 
interaction  potentials,  the  surface  coupling  matrix  elements,  and  the 
amount  of  depletion  from  competing  continuum  coupling  and  excitation 
transfer  channels.  With  only  crude  model  potentials  and  order-of-magnitude 
estimates  for  the  continuum  coupling  and  excitation  transfer  effects,  it 
does  not  seem  justified  to  vary  parameters  to  achieve  a  fit  to  the  data 
that  would  be  both  model -dependent  and  non-unique.  However,  an  estimate 
of  the  probable  causes  for  the  differences  from  reaction  (1)  can  be  useful 
in  suggesting  directions  for  further  experimental  or  theoretical  effort 


on  this  system. 


The  shift  of  the  ionic  rainbow  for  reaction  (1),  to  lower  T  values 

than  that  for  reaction  (2),  seems  to  be  consistent  with  a  shallower  well  in 

the  Ar+  +  interaction  (see  Figure  6).  This  change  also  lowers  the 

magnitude  of  the  covalent  peak  relative  to  the  ionic  peak.  In  fact  all 

other  modifications  combined  would  only  need  to  lower  the  covalent  peak 

by  an  additional  30-40%  to  achieve  a  good  agreement  with  the  experimental 

data  at  E  “  51  eV.^ 

0 

Combined  with  the  change  in  the  ionic  potential  described  in  Section 

E,  an  increase  of  by  -407.  over  the  K  +  value  (0.28  eV)  would  achieve 

a  reasonable  match  to  the  experimental  data.  However,  such  a  change  in  V 

still  seems  rather  large.  As  evidence  for  this  assertion,  we  note  that 
26 

Hubers  et  al.  present  a  reduced  matrix  element  formula  that  is  slightly 

modified  from  Eq.  (7)  and  correlates  with  all  alkali-halogen  systems  to 

* 

an  accuracy  of  better  chan  10%.  The  modifi  ation  of  the  Ar  +  matrix 

26 

element  postulated  here  would  deviate  by  40%  from  the  prediction  based 

on  the  alkali  data  and  would  then  suggest  some  real  differences  of  Ar* 

from  alkali-like  properties  at  distances  of  R  «  5.7  A.  It  is  difficult 

c 

Co  understand  such  a  large  effect  at  that  distance  caused  by  a  vacancy 
in  the  ion  core. 

Alternatively,  a  continuum  coupling  width  of  —40  T  could  be  combined 

o 

with  the  shallower  ionic  potential  to  match  the  peak  ratios  adequately. 

This  width  is  somewhat  larger  chan  the  rough  estimate  (10-15  T  )  based 

0 

10 

on  the  Miller-Morgner  arguments,  and  the  attenuated  covalent  peak  would 


be  located  at  larger  t  values  than  the  data  would  allow.  Yet  our  estimate 
of  T  is  subject  to  considerable  uncertainty  and  could  be  significantly 
in  error.  Also,  the  covalent  peak  position  and  shape  are  sensitively 

4 

dependent  on  the  potential  functions  used.  This  has  been  demonstrated 
in  K  +  Br^  by  a  comparison  of  calculations  on  diabatlc  and  adiabatic 
surfaces.  Hence,  continuum  coupling  (3)  and  excitation  transfer  (27) 
reactions  involving  the  ten  valence  electrons  could  be  responsible  for 
all  or  part  of  the  remaining  30-403  depletion  of  the  covalent  peak. 

Finally,  we  note  that  making  the  potentials  more  repulsive  tends  to 
decrease  the  size  of  the  covalent  peak  relative  to  the  ionic  peak.  For 
simplicity,  we  have  constrained  our  calculations  so  that  the  repulsive 
Born-Mayer  terms  are  identical  on  the  ionic  and  covalent  surfaces. 

Although  this  is  a  reasonable  first  approximation,  one  could  decrease  the 
covalent  peak  by  relaxing  this  constraint  and  allowing  more  repulsion 
in  the  covalent  potential.  If  such  a  difference  were  responsible  for 
the  smaller  covalent  peak  in  (1)  relative  to  (2),  it  would  necessarily 
imply  differences  between  the  two  systems  in  the  repulsive  region. 

This  is  easier  to  accept  than  large  differences  (403)  in  the  coupling 
matrix  elements  at  the  much  larger  distances  associated  with  the 
surface  crossings. 
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IV  SUMMARY 


The  small  differences  between  the  differential  cross  sections  for 

* 

ion-pair  formation  in  the  Ar  +  and  R  +  1^  system  can  be  explained, 

^  • 

in  part,  by  a  more  repulsive  wall  in  die  ionic  Ar  +  surface  relative 

to  K  +  I2  .  This  yields  a  smaller  ionic  well  depth  (by  -0.5  eV  at  r^ 

of  I^)  and  a  correspondingly  smaller  ionic  rainbow  angle  that  is  in  agree- 

* 

mant  with  the  data.  If  the  covalent  surfaces  Ar  +  I,  and  K  +  I  have 

2  2 

repulsive  terms  similar  to  those  on  their  ionic  counterparts,  the  calculated 

ratio  for  the  intensity  of  the  covalent  peak  to  the  ionic  peak  in 

reaction  (1)  is  too  large  by  a  factor  of  ~1.4. 

This  ratio  can  be  improved  by  one  or  more  of  the  following:  a 

* 

relative  increase  in  the  repulsion  on  the  Ar  +  covalent  surface;  a 
depletion  of  the  covalent  peak  by  competing  Penning  (3)  and  excitation 
transfer  (27)  reactions;  and  an  increase  in  the  coupling  matrix  element 
at  the  surface  crossing  between  the  ionic  and  covalent  surfaces.  The 
third  possibility  seems  least  likely  since  it  requires  effects  from  the 
Ar  core  vacancy  to  be  important  at  large  R  (—5.7  X) . 

Since  the  potential  surfaces  have  only  been  estimated  in  a  crude 
way  and  since  the  continuum  coupling  width  is  only  known  very  approxi¬ 
mately,  we  feel  that  no  claim  to  a  unique  fit  can  be  made.  Continuum 
coupling  and  excitation  transfer  processes  may  Indeed  cause  observable, 
but  small,  depletions  of  the  covalent  peak  in  the  ion-pair  process  (1). 
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FIGURE  CAPTIONS 


1.  Two  possible  paths  a  and  b  to  ion  pair  production  are  schematically 

* 

indicated.  The  potential  curves  represent  cuts  of  the  Ar  +  I  and 

+  2 
Ar  +  I^  surfaces  at  a  fixed  I^  separation. 

2.  Illustration  of  the  coordinate  system  used  to  define  the  potential 
functions. 

3.  Comparison  of  previous  calculations  and  experimental  data  (Ref.  2) 
on  the  system  K  +  with  the  present  calculation,  which  uses 
slightly  different  interaction  potentials. 

4.  Calculation  of  the  differential  cross  section  for  ion  pair  formation 

* 

in  Ar  +  I^  compared  with  the  data  of  Ref.  1.  The  effect  of  different 
autoionization  rates  T(R)  in  the  covalent  channel  is  illustrated, 
r  is  the  function  estimated  using  the  method  discussed  in  Section 
II-B-1,  and  given  by  Eq.  (26).  Successively  larger  values  of  I"(R) 
reduce  the  covalent  peak  by  larger  fractions. 

3.  Calculation  of  the  differential  cross  section  for  ion  pair  formation 
* 

in  Ar  +  I 2  compared  with  the  data  of  Ref.  1.  The  effect  of  different 
coupling  matrix  elements  V ^  Is  illustrated.  T(R)  is  set  to  zero  in 

all  cases.  V°  is  the  form  given  by  Eq.  (7)  with  c  -  0.28  eV  and 

•-1  1 
c^  «  0.65  A  ,  the  same  parameters  used  in  Ref.  1  for  K  +  I^.  The 

scaled  values  of  correspond  to  changing  c  to  0.7  eV  and  1.4  eV. 

6.  Calculation  of  the  differential  cross  section  for  ion  pair  formation 
* 

in  Ar  +  I^  compared  with  the  data  of  Ref.  1.  The  potential  parameters, 
defined  in  Eqns.  (4)  -  (6),  are  A  »  6.4  x  104  eV  and  B  ■  4.762  A"1  for 
the  original  well,  and  A  -  1.35  x  104  eV  and  B  -  3.781  A"1  for  the 
shallower  well. 
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MODEL  FOR  FAST,  NON -ADIABATIC  COLLISIONS  BETWEEN 
ALKALI  ATOMS  AND  DIATOMIC  MOLECULES 

A.  P.  Hickman 

Molecular  Physics  Laboratory 
SRI  International 
Menlo  Park,  CA  94025 

ABSTRACT 

Equations  for  collisions  involving  two  potential  surfaces  are 
presented  in  the  impact  parameter  approximation.  In  this  approxi¬ 
mation,  a  rectilinear  classical  trajectory  is  assumed  for  the 
translational  motion,  leading  to  a  time -dependent  Schroedinger ' s 
equation  for  the  remaining  degrees  of  freedom.  Model  potentials  are 
considered  for  collisions  of  alkali  atoms  with  diatomic  molecules 
that  lead  to  a  particularly  simple  form  of  the  final  equations. 

Using  the  Magnus  approximation,  these  equations  are  solved  for  parameters 
chosen  to  model  the  process  Cs  +  0^  -  Cs+  +  0^  ,  and  total  cross  sections 
for  ion-pair  formation  are  obtained  as  a  function  of  energy.  The  results 
exhibit  oscillations  that  correspond  qualitatively  to  those  seen  in 
recent  measurements.  In  addition,  the  model  predicts  that  the  oscilla¬ 
tions  will  become  less  pronounced  as  the  initial  vibrational  level  of 
0^  is  Increased. 
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I .  INTRODUCTION 


In  ehis  paper  we  consider  Che  process  of  ion-pair  formation  in 

collisions  between  alkali  atoms  and  diatomic  molecules: 

A  +  XY  -  A+  +  XT’  (1) 

An  important  feature  of  this  process  is  that  two  electronic  potential 

surfaces  are  involved,  the  location  of  whose  crossing  depends  both  on 

the  vibrational  and  translational  coordinates.  Considerable  work  has 

been  done  to  analyze  experimental  data  in  terms  of  a  perturbative, 

1-3 

classical  trajectory,  surface-hopping  model.  The  agreement  is 
satisfactory  for  many  systems.  However,  there  has  been  a  tendency  to 
try  to  improve  the  basic  model  by  introducing  various  ad  hoc  modifica¬ 
tions.  Here  we  will  attempt  to  treat  the  process  using  a  more 
systematic  theory. 

Our  approach  is  based  on  the  semiclassical  impact  parameter 
4 

approximation.  We  choose  a  model  form  for  the  two  potential  surfaces 
that  is  consistent  with  previous  work,^’^  and  also  leads  to  a  particularly 
simple  form  of  the  final  equations  of  motion.  We  assume  that  the 
translational  motion  is  a  straight  line  trajectory  with  constant  velocity, 
and  then  focus  on  the  quantum  mechanical  evolution  of  the  remaining 
degrees  of  freedom  under  the  resulting  time-dependent  Hamiltonian.  An 
approximate,  unitary  S  matrix  is  found  using  the  Magnus  approximation.^ 


In  Section  II  we  present  the  equations  of  notion  in  the  impact 


parameter  approximation  for  collisions  involving  two  potential  surfaces 
that  are  functions  of  vibrational  and  translational  coordinates.  A 
matrix  notion  is  used  chat  highlights  the  manner  in  which  the  second, 
potential  surface  affects  the  motion.  It  is  shown  how  rotational 
degrees  of  freedom  may  be  created  using  the  sudden  approximation.  In 
Section  III  results  are  obtained  for  potential  surfaces  chosen  to  model 
the  specific  system 

Cs  +  02  -  Cs+  +  0  \  (2) 

Total  cross  sections  for  ion-pair  formation  as  a  function  of  energy 

are  calculated.  Because  of  the  model  nature  of  the  potential  surface, 

and  the  perturbation  solution  of  the  equations  of  motion,  the  results 

are  not  expected  to  be  quantitative.  However,  the  results  are  in 

8 

qualitative  agreement  with  recent  experiments.  Concluding  remarks 
are  given  in  Section  IV. 


\ 


h 
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II.  THEORY 


A.  General  Equations 

We  first  develop  the  equations  of  the  impact  parameter  approximation 
for  the  case  of  tvo  electronic  potential  surfaces  that  depend  on  a 
translational  coordinate  R  and  a  vibrational  coordinate  r.  In  the  next 
section  we  shall  derive  Che  form  taken  by  these  equations  for  our 
particular  model  potential,  and  show  how  the  rotational  motion  may  be 
included  using  the  sudden  approximation. 

Let  us  cake  the  matrix  elements  of  the  total  Hamiltonian  of  the 
system  in  a  basis  of  two  (diabetic)  electronic  functions.  We  will 
assume  chat  che  resulting  operator  has  the  following  form: 


H 


(Tr+T_).1  + 


Vr>+ VR) 


V1(r)  +  Vl(R). 


VR’r)  VR’r) 

W01(R,r)  wu(R.r> 


(3) 


Mote  chat  we  have  isolated  the  separable  part  of  che  diagonal  potential 
matrix  elements.  The  translational  and  vibrational  coordinates  are  R 
and  r,  respectively.  We  now  assume  chat  che  translational  motion  is  given 
by 

fz  2> 

R(C)  -  b  +(vt)4  (4) 

where  b  and  v  are  che  impact  parameter  and  relative  velocity,  respectively. 


We  can  now  write  che  following  time -dependent  Hamiltonian  for  the 


(t)  -^(r)|W0i(R(t),ir)|4^(r)^  exp  i.  J  £e°  +  ^(tS-e^-a^Ct ' )J  dt ' 


A  is  the  Hemitian  conjugate  o£  A: 


It  , 

-i  J  £  e°  +  -JjQ(t •  )-«^-o>1(t '  )J  dt'j 


There  will  be  2N  solution  vectors,  corresponding  to  the  system  beginning 

at  t  •  -®  in  each  of  the  possible  states  j  ,  y.»l,2 . N  and  j$^j, 

X  *  l...,N.  If  we  solve  Eq.  (12)  subject  to  the  initial  condition  that 
the  solution  matrix  is  the  unit  matrix  at  t  -*  -®  ,  then  the  value  of  the 
solution  at  t-*+®  can  be  used  to  find  the  probability  that  the  system 
is  in  each  possible  final  state  after  the  collision.  Each  element  of 
the  solution  matrix  at  t  ■  +*  is,  within  a  phase  factor,  equal  to  the 
corresponding  element  of  the  S  matrix.  This  phase  factor  will  be 
discussed  later. 


B.  Application  to  a  Model  System 
1.  Isotropic  Potentials 

In  this  section  we  will  consider  the  form  of  Eq.  (12)  for  a 
particular  set  of  potentials  chosen  to  model  collisions  of  alkali  atoms 
with  diatomic  molecules.  Such  systems  are  characterized  by  surface 
crossings  at  fairly  large  values  of  R  (>6a^  ),  so  that  near  the  crossings 
fairly  simple  analytic  approximations  may  be  appropriate.  We  will  use 
analytic  potentials  similar  to  those  already  discussed  in  previous  work:^' 
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V0(R)  -  0  (17) 

V  (R)  -  AE  -  —  (18) 

1  R 

woo(R)  “  wu(R)  “  0  (19) 

W01(R>  *  C1  exP(c2R)  •  (20) 

vt(r)  -  aJ  exp  £  -2@t(r  -  r)  J  -2  exp^-g^r  -  r^j!  (21) 


These  analytic  functions  model  the  behavior  of  the  surface  at 
large  values  of  R.  The  covalent  surface  is  approximately  constant,  the 
ionic  surface  is  coulombic,  and  their  asymptotic  separation  is  AE.  Such 
potentials  are  expecced  to  be  adequate  as  long  as  the  major  part  of  the 
cross  section  comes  from  trajectories  with  large  impact  parameters.  It 
should  be  noted  Chat  one  could  add  the  same  repulsive  term  to  VQ  and 
without  affecting  the  transition  probabilities,  since  the  formulas  to  be 
presented  here  depend  only  on  Che  difference  potential.  This  modification 
would  allow  Vq  and  to  be  more  realistic  at  small  R. 

Equations  (12)-(21)  lead  to  Che  following  simplifications  in  Eq.  (12). 

t 

First,  the  diagonal  blocks  F  and  G  vanish  entirely  because  of  Eq.  (19). 

'  These  are  Che  terms  that  lead  to  vibracionally  Inelastic  scattering  on  a 

i 

4  single  surface.  Second,  the  matrix  elements  of  W^^(R(t),r)  reduce  Co  a 

function  of  c  times  the  Franck-Condon  factor 

CuX  CM) 


U 

i 

'¥ 
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because,  according  Co  Eq,  (20),  W  is  independent  of  r.  Finally,  Che 
integral  inside  the  exponential  can  be  done  analytically.  The  final 
form  of  the  equations  is  thus 


where 


V<c) '  Voi(c)  e!tp  llf(c)i 

(24) 

*01(C)  ’  C1  exp  |-c2[b2+ (Vt)2j  j- 

?25) 

f(t)  -  -(«+«£  -«“)c+iln|f  +[l+(f)] 

1 1  (26> 

In  many  applications  it  will  probably  be  necessary  to  solve  Eq.  (23) 
numerically.  However,  in  the  present  paper  we  will  resort  to  an 


approximation  that  involves  considerably  less  computation.  We  use  the 
Magnus  approximation,2  which  involves  setting  the  solution  U  at  t  -*  +• 
to  the  exponential  of  a  matrix.  The  result  is 
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where  U  is  a  2N  x  2N  matrix.  The  matrix  Q  and  its  transpose  Q  are  defined  by 


09 

L 


V*  /  Vc,d' 


2C 


aa 

J  ju01(t)  C0s[f(c 


(28) 


)]  dt 


Because  o€  the  form  of  Eq.  (27),  unitarity  is  preserved.  In  the  limit  of 

weak  coupling,  one  may  expand  the  exponential  of  Eq.  (27)  and  recover  the 

standard  result  of  first  order  time  dependent  perturbation  theory.  Note 

also  that  integral  in  Eq.  (28)  could  be  evaluated  using  the  stationary 

phase  approximation.  The  integrand  is  rapidly  oscillating  except  near  the 

time  corresponding  to  the  surface  crossing,  so  that  the  formalism  reflects 

the  expected  semiclassical  result.  In  the  present  work,  however,  the 

integrals  will  be  evaluated  by  numerical  quadrature. 

We  will  denote  the  elements  of  U  by  (b),  where  0  £  i,j  £  1 

and  1  £  X.p.  £  N.  The  superscripts  denote  the  NxN  blocks,  the  subscripts 

the  elements  within  each  block.  In  this  paper,  we  are  primarily  interested 
10 

in  U.  (b),  the  block  of  transition  amplitudes  between  an  initial  covalent 

A.H 

channel  and  a  final  ionic  shannel.  For  this  case,  the  total  cross  section 


from  an  initial  vibrational  state  p,  to  a  final  state  \  is  given  by 


We  may  also  define  a  cross  section  summed  over  final  states: 

N 

ff(u>  -  .1,  9  (30) 

1 

We  now  briefly  discuss  the  determination  of  the  S  matrix  from  the 

matrix  U.  Recall  that  the  squared  magnitude  of  each  element  of  U  gives 

the  probability  of  the  corresponding  state-to-state  transition.  That  is, 

each  element  of  U  must  be  multiplied  by  some  phase  factor  in  order  to 

obtain  S.  This  phase  factor  is  not  necessary  for  the  present  applications 

[ Eqs .  (29)  and  (30)],  but  it  would  be  needed  for  a  differential  cross 

9 

section.  We  adopt  the  procedure  of  multiplying  elements  of  U  by  phase 
factors  depending  on  ^  and  T]  ,  which  are  the  single  channel  phase  shifts 
obtained  for  the  potentials  Vq(R)  and  V^(R)^  respectively.  [Note  that 
7]  is  the  phase  shift  in  addition  to  the  standard  coulomb  phase  shift] . 

The  result  may  be  summarized  by  writing 

S  -  A  0  A  (31) 

*>*«*  *■*«*  aw 


where  A  is  a  2N  x  2N  matrix  given  by 


(32) 


This  formal  prescription  amounts  to  multiplying  each  element  of  the  ij 
block  of  £  by  exp  ^1(7^  +  H ^)J  ,  fori,J-0,l. 
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2.  Anisotropic  Potentials 


To  include  rotational  degrees  of  freedom  in  the  model  it  is  necessary 
to  include  the  anisotropy  of  the  interaction  potential.  We  generalize  the 
potentials  in  the  following  way,  according  to  the  coordinates  shown  in 
Figure  1. 

Vq(R,Y)  -  V(Rl)  +  V(R2)  (33) 

V,(R,y)  -  V(R,v)  +  AE  -  ~  “  (34) 

1  Rl  *2 

1  a2 

-  V  (R,y)  +  AE  -  Pfcos  v)  +  ...  (35) 

0  R  R3  2 

"C2R 

W  (R*y)  *  c:  e  *  sin  2y  (36) 

The  potential  V  is  an  effective  two-body  potential  between  the  projectile 

and  each  atom  of  the  (homonuclear)  diatomic  target,  which  is  assumed  to 

have  an  intemuclear  separation  of  2a.  As  in  the  isotropic  case,  an 

exponential,  repulsive  form  of  V(R)  could  be  chosen.  However,  for  the 

particular  form  of  Eqs.  (33)  and  (34),  the  results  to  be  presented  here  do 

not  depend  on  V  since  it  does  not  appear  in  the  difference  potential. 

This  form  of  and  has  been  used  previously.^’^  The  angular  dependence 

8 

of  W^  is  chosen  to  be  the  same  used  by  Kleyn  et  al.  for  Cs  +  0^,  and  has 

10 

nodes  at  the  angles  required  by  symmetry. 

We  now  obtain  rotationally-averaged  cross  sections  using  the  Inflnlte- 

11  12 

Order-Sudden  (IOS)  approximation,  ’  in  which  the  cross  section  averaged 


D-12 


over  initial  rotational  states  and  summed  over  final  states  is  obtained  by 
averaging  the  cross  sections  obtained  in  a  sec  of  calculations  in  which  y 
is  regarded  as  a  parameter.  In  ocher  words,  the  rotacionally  averaged 
cross  section  from  an  initial  vibrational  state  to  a  final  vibrational 
state  \  is  given  by 

1 

%  "  2 

0 

where  a.  (y)  is  the  cross  section  obtained  for  "frozen"  y.  When  y  is 
fixed  ic  is  easy  to  calculate  the  cross  sections  for  ion  pair  formation, 
using  the  methods  of  the  preceeding  section,  for  the  form  of  VQ  and 
given  in  Eqs.  (33)  and  (34).  The  result  for  0(y)  is  given  by  an  equation 
of  Che  same  form  as  Eq.  (27),  with 


TT 

/ 


o.  (Y)  sin  ydy 


(37) 


where 


Q  (y)  -  2C  sin2y 

ILA. 


OB 

j  «01Ct)  "*[fY(e)] 


dt 


a  P  (cosy) 
f^c)  -  f(t) - - 


vb 


(38) 


(39) 
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The  formalism  developed  in  che  preceeding  section  has  been  applied 
with  parameters  chosen  to  model  the  reaction 

Cs  +  02  -  Cs+  +  0*  .  (40) 

Initially  we  used  an  interaction  potential  without  anisotropy.  In  a 
second  calculation,  however,  anisotropy  is  included,  and  rotationally 
averaged  cross  sections  are  obtained  using  the  IOS  approximation.  For 
both  of  the  calculations  we  included  12  vibrational  states  of  0^  and  12 
of  O'.  This  number  was  sufficient  for  convergence.  The  energy  levels 

e*  and  and  overlap  integrals  C  were  calculated  numerically  using 

U  \  ’M 

the  Morse  potentials  for  0^  and  0^  defined  by  Eq.  (21)  and  Table  1. 

Other  parameters  for  the  potentials  are  al3o  listed  in  Table  1. 

The  results  of  the  first  calculation  are  shown  in  Figure  2.  Total 

cross  sections  for  ion-pair  formation  are  shown  as  a  function  of  energy. 

The  cross  sections  are  summed  over  transitions  from  the  indicated  initial 

vibrational  state  to  all  final  vibrational  states.  The  most  interesting 

feature  is  che  pronounced  oscillations.  These  have  been  previously 

2,6 

interpreted  in  terms  of  a  classical  trajectory-hopping  model,  which 
is  briefly  summarized  here.  The  basic  idea  is  chat  the  crossing  point 
(or  seam)  of  the  covalent  and  ionic  potential  surfaces  depends  on  the 


vibrational  coordinate.  On  the  inward  part  of  the  trajectory,  if  the  0^ 
is  in  the  v»0  state,  the  vibrational  coordinate  is  localized  near  its 
equilibrium  value,  and  the  crossing  occurs  at  Ra7.2  a^,  If  a  crossing 
to  the  ionic  surface  occurs,  however,  0^  will  be  formed,  and  the 
vibrational  coordinate  will  start  to  expand,  since  the  equilibrium 
separation  of  0^  is  larger  than  that  of  0^.  Depending  on  the  relative 
translational  velocity,  the  0^  may  undergo  a  fraction  of  a  vibration  or 
several  vibrations  before  the  crossing  seam  is  reached  a  second  time. 

If  the  0^  is  fully  expanded  at  the  second  crossing,  this  crossing  will 
occur  at  a  larger  value  of  R,  and  the  matrix  element  for  recrossing  to 
the  covalent  surface  will  be  much  smaller,  because  of  the  exponential 
dependence  of  WQ^(R)  •  In  other  words,  if  the  02  undergoes  a  half-integral 
number  of  vibrations  between  the  surface  crossings,  the  cross  section  for 
ion-pair  formation  will  be  enhanced. 

This  classical  model  has  successfully  explained  the  oscillations 

8 

observed  in  the  experimental  data  of  Kleyn  et  al.  It  is  worth  noting, 

13 

however,  that  except  for  a  few  exceptions,  the  previous  calculations 
have  assumed  that  the  0^  begins  at  rest  in  the  bottom  of  its  potential 
well.  Our  present  calculations  show  chat  the  oscillations  persist  in  a 
model  that  treats  the  diatomic  target  quantum  mechanically.  He  are  also 
able  to  examine  the  effect  of  changing  the  initial  vibrational  state  of 
02>  It  is  seen  that  the  oscillations  tend  to  become  less  pronounced  as 
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v^  is  increased.  A  possible  explanation  for  this  effect  is  that  as  v^ 
increases,  the  spatial  extent  of  the  0^  vibration  increases,  and  it 
becomes  more  difficult  to  localize  the  first  crossing.  The  oscillations 
depend  on  the  first  and  second  crossings  occurring  at  very  different 
values  of  R. 

In  Figure  3,  we  show  the  results  of  calculations  in  which  an 
anisotropic  term  of  the  potential  is  treated  using  the  IOS  approxima¬ 
tion.  Compared  to  Figure  2,  it  is  seen  that  the  anisotropy  causes 
the  oscillations  to  be  less  pronounced,  but  that  the  same  trends  persist. 
Figure  4  compares  the  theoretical  results  for  v  ■  0  with  the  data  of 
Kleyn  et  al.  The  trajectory  hopping  calculation  of  Kleyn  et  al.6  is 
also  shown.  All  theoretical  cross  sections  have  been  multiplied  by  a 
'statistical  factor  of  1/3.  Qualitative  agreement  between  experiment 
and  the  present  theory  is  obtained.  It  is  likely  that  the  position  of 
the  maximum  in  the  theory  could  be  brought  into  closer  agreement  with 
experiment  by  adjusting  the  parameters  c^  and  c^  of  the  coupling.  He 
have  simply  taken  the-  values  obtained  by  Kleyn  et  al.6  by  fitting  their 
classical  calculations  to  the  data.  (If  the  data  of  Kleyn  et  al. 
presented  in  Figure  4  are  scaled  by  a  constant  factor,  the  agreement 
with  the  classical  calculation  is  reasonable.  This  procedure  is  justi¬ 
fied  because  the  experimental  cross  sections  are  relative.) 
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A  possible  refinement  to  the  theory,  which  would  improve  the 

threshold  behavior,  would  be  to  include  a  factor  (E^/E^)  1°  the  cross 

section  formulas.  An  empirical  correction  of  this  form  has  sometimes 
14 

been  used  to  insure  that  the  cross  section  goes  to  zero  at  threshold. 
For  energies  near  threshold,  the  assumption  of  a  straight-line,  constant 
velocity  trajectory  is  not  valid.  We  have  not  included  this  factor  in 
the  present  calculations.  For  the  cross  sections  at  the  lowest  energies 
presented  in  Figures  2,  3,  and  4,  the  correction  would  be  on  the  order 
of  10%. 

It  is  Interesting  to  consider  what  might  cause  the  difference 
between  the  surface-hopping  calculation  and  the  present  model.  Both 
use  essentially  the  same  potentials,  yet  the  difference  at  some  energies 
is  as  much  as  a  factor  of  two.  The  discrepancy  may  be  related  to  the 
following  points.  First,  the  rotation  of  the  target  is  handled 
differently  in  the  two  models.  In  the  present  work,  introducing  the 
IOS  approximation  causes  the  location  of  the  broad  peak  in  the  cross 
section  to  shift  about  40%  (cf.  Figures  2  and  3).  This  significant 
change  suggests  the  cross  sections  may  be  sensitive  to  the  treatment  of 
the  molecular  anisotropy.  Second,  the  Magnus  approximation  may  not 
provide  a  sufficiently  accurate  solution  to  Eq.  (23).  Finally,  the 
Landau-Zener  formula  used  in  the  classical  model  often  underestisiaces 
transition  probabilities  at  high  energies.1^  Further  work  is  necessary 
before  firm  conclusions  can  be  drawn. 


T  \ 
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IV.  CONCLUDING  REMARKS 


We  have  presented  Che  equations  of  the  semiclassical  impact 
parameter  approximation  for  the  case  of  collisions  involving  two 
potential  surfaces,  and  have  examined  in  detail  their  application  to 
collisions  of  alkali  atoms  with  diatomic  molecules.  Our  results  for 
the  specific  system  Cs  +  0^  are  encouraging  and  suggest  that  the  method 
may  be  of  great  value  in  fast  collisions  where  successive  curve  crossings 
are  not  well  isolated. 

The  present  application  has  depended  heavily  on  rather  specific 
assumptions  about  the  analytic  form  of  the  interaction  potential.  These 
assumptions  were  useful  in  order  to  simplify  the  numerical  work  required 
in  the  first  application  of  the  theory.  However,  a  number  of  refinements 
may  be  introduced  in  a  manner  that  is  still  computationally  feasible. 

In  particular,  more  general  forms  of  the  potentials  V^(R)  and  V^(R)  may 
easily  be  Introduced.  Including  the  non-separable  terms,  Vqq(R>7)  and 

g 

W^(R,r)  is  somewhat  more  difficult.  Fortunately,  experimental  evidence 
seems  to  suggest  that  these  terms  may  not  be  too  important  for  collisions 
of  alkali  atoms  and  diatomic  molecules. 

Future  work  will  focus  on  the  applicability  of  the  10S  to  collisions 
involving  two  potential  surfaces,  and  on  the  application  of  more  refined 
numerical  techniques  for  solving  the  equations  of  motion. 
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Figure  Captions 


Figure  1.  Diagram  illustrating  the  coordinates  used  on  the  present 
calculations. 

Figure  2.  Results  o£  the  model  calculation  of  the  cross  section  for 
ion-pair  formation  in  collisions  of  Cs  +  0^.  An  isotropic  potential 
is  used  as  defined  in  Eqs.  (17)-(21)  and  Table  1.  The  cross  sections 
are  summed  over  all  (12)  final  vibrational  states. 

Figure  3.  Results  of  the  model  calculation  of  the  cross  section  for 
ion-pair  formation  in  collisions  of  Cs  +  0^.  The  anisotropic  potential 
defined  by  Eqs.  (33)-(36)  is  used,  and  the  cross  sections  are  summed 
over  all  (12)  final  vibrational  states. 

Figure  4.  Comparison  of  the  present  model  with  the  experimental 
results  and  calculations  of  Kleyn  et  al.  The  latter  calculation  was 
fitted  to  the  data;  if  the  experimental  cross  sections,  which  are 
relative,  are  scaled  by  a  constant  the  agreement  with  that  calcula¬ 
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